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ABSTRACT 


This dissertation investigates power requirements for a Free Electron Laser to 
bum through various missile radome materials. It also includes computer simulation 
results for several FEL system configurations designed to achieve maximum power while 
maintaining strict energy spread constraints. 

The method used to determine power requirements to bum through materials was 
to use the Thomas Jefferson National Accelerator Facility’s Free Electron Laser to 
conduct material damage experiments. As the laser was improved and increased in 
power, the laser spot sizes on the target materials were increased while maintaining a 
constant irradiance. The key results from these experiments included determining 
minimal spot sizes that can be used for future experiments, and validation that an 
irradiance level of 10 kW/cm^ can bum through most missile radome materials in a few 
seconds. 

The computer simulations involved changing various parameters of a FEL such as 
electron energy levels, pulse lengths, magnetic field strengths, desynchronism, as well as 
several other parameters, to determine the best possible configuration to achieve the 
desire power levels and energy spread requirements for development of a megawatt size 
FEL. The results indicate that for the proposed designs, both the required power and the 
required energy spread limit can be met. 
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I. 


INTRODUCTION 


Since the end of the Cold War the United States Navy 
has adopted operational strategies that follow the ideas 
presented in From the Sea and Forward. . .From the Sea [Refs. 
1 and 2] . The central idea in those dociiments is the 
necessity of a new maritime concept, one in which the Navy 
plays a fundamental role in shore operations being carried 
out by our airmed forces. As a result Navy ships, 
submarines, and aircraft can be expected to operate in the 
littoral regions of enemy territories, with the increased 
potential for encoxmtering hostile forces. Even with the 
advanced defensive systems on board most of our surface 
ships, the decreased reaction times available in littoral 
operations makes them considerably more vulnerable to anti- 
ship cruise missiles (ASCMs) and attacks by fast attack 
boats. Additionally, the close proximity of hostile, 
neutral, and friendly forces makes precision targeting 
extremely important. 

Building on this new maritime concept, the Chief of 
Naval Operations recently released the Navy Strategic 
Planning Guidance (April 2000) which prioritizes the 
capabilities needed to support such a maritime concept. The 
priority levels range from Priority (I), "Those capabilities 
that directly support or enhance the enduring core naval 
competencies without which severe strategic risk would be 
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incurred", down to Priority (V), "Those capabilities that 
directly support or enhance... without which minimal 
strategic risk would be incurred." [Ref. 3, pgs. 40, 61 and 
62] 

Listed under priority (I) capabilities are: 

• The Fleet-wide, point defense capability to achieve 
high-probability hardkill against sub- and super¬ 
sonic cruise missiles. 

• The capability to provide area air and missile 
defense against emerging threats, including advanced 
cruise missiles. 

• The capability to detect, identify, track and 
destroy high numbers of small craft in the 
littorals. 


According to the Office of Naval Intelligence (ONI), 
the trends in ASCMs will include: 

• Significant increase in missile terminal velocity 
from predominately sub-sonic speeds today to 
predominately supersonic by 2010 and hypersonic by 
2020. 

• Complex terminal maneuvers and seekers to tax point 
defense systems. 


Additionally, ONI concludes that the capabilities most 
likely to be encountered include massed small boat attack, 
unalerted single salvo cruise missile attack, short-range 
ballistic missiles, and terrorist attacks. [Ref. 3, pgs. 8 
and 14] Nearly 70 nations have deployed sea and land 
launched ASCMs, and 20 nations have air launched ASCMs [Ref. 
4, pg. 4]. 
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Super-sonic and highly maneuverable missiles are 
already held by known hostile nations, including the Russian 
built SSN-22 Sunburn which has a cruise speed/terminal speed 
of mach 2.5/mach 4.5 [Ref. 5, pg. 121]. Though the terminal 
speed alone is probably enough to greatly reduce the 
effectiveness of most self defense systems, this missile 
also includes a terminal maneuver. This will reduce the 
effectiveness of self-defense systems such as the Navy's 
Close in Weapon System (CIWS), which is especially 
vulnerable to these maneuvers [Ref. 6, pg. 179]. 

One of the nations that hold this missile is Iran, 
which has the land launched version [Ref. 5, pg. 325] . If 
laiinched in the Straits of Hormuz, which is approximately 40 
nautical miles wide, it would take less than 50 seconds to 
impact a ship transiting approximately 20 nautical miles 
away through the central section of the straits. Current 
self-defense systems would have an extremely difficult time 
engaging such a threat within the 50 second window, 
especially if it was an unalerted launch. It is threats 
such as these that make the idea of a Directed Energy Weapon 
(DEW) stand out as a possible solution. 

One such DEW would be a High Energy Laser. In a study 
carried out by the National Research Council to determine 
Technology needed for the Navy and Marine Corps for the 
years 2000 - 2035, a tunable Free Electron Laser (PEL) is 
specifically called out as the laser of choice for such a 
weapon due to its ability to be tuned to wavelengths that 
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offer minimal atmospheric absorption and scattering [Ref. 6, 
pg. 184] . A secondary advantage of using a FEL is that its 
gain medium is a vacu\im, therefore there is no requirement 
for the removal of large amounts of waste heat from the 
lasing mediiim as with a solid state laser, or the removal of 
toxic effluents associated with chemical lasers. 

Several topics are presented in this dissertation on 
the use of a FEL as a Directed Energy Weapon. Chapter II 
discusses the results of material damage experiments 
conducted at Thomas Jefferson National Accelerator Facility 
(TJNAF) using a FEL. Some of this work was originally 
presented by the author at the 21®’' International FEL 
Conference in Hamburg, Germany, and is expected to be 
piiblished in Ref . 7 . 

Chapter III discusses classical FEL Theory. Starting 
with the Lorentz force equation and Maxwell's wave 
equation, the dimensionless FEL pendulum and wave equations 
are developed [Ref. 29]. The work in this thesis includes 
modifications of the FEL interactions and gain to include 
short pulse effects and tapering are also developed, 
including a modified pendulum equation to include the 
effects of tapering of the undulator. This work was 
presented by the author at the 21®" International FEL 
Conference in Hamburg, Germany and was published in Ref. 8. 
More recent results from simulations conducted by the author 
were presented at the 22®^ International FEL Conference in 
Raleigh, NC, and will be published in Ref. 7. 
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Chapter IV includes the results of computer simulations 
for proposed output power upgrades to the TJNAF FEL. These 
upgrades are the next step in a proposed upgrade to 1 MW 
output power. This work was analyzed by the author and was 
presented at the 20"*' International FEL Conference in 
Williamsburg, VA, and 22"^ International FEL Conference in 
Raleigh, NC. Results from 20 kW simulations were published 
Ref. . 9, and results for the 10 kW simulations will be 
published in Ref. 7. 

The last chapter discusses the effect of tapering the 
undulator, and specifically calls out the differences and 
advantages of positive vs. negative taper. Simulations 
conducted by the author resulted in surprising positive 
results, and were presented by him at the 21*‘ and 22“^ 
International Free Electron Laser Conferences, and published 
in Refs. 7 and 8. 
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II. MATERIAL DAMAGE EXPERIMENTS 


A. THEORY 

The interaction between lasers and matter is a 
complicated issue generally described by non-linear and 
unpredictable effects. The basic kill mechanism for 
destroying a missile is to cause a catastrophic failure of 
one of the central components such as the warhead, control 
and navigation electronics, or the high pressure chamber in 
the propulsion section, or to create aerodynamic instability 
beyond the design limits of the missiles flight control 
surfaces. [Ref. 10] 

The laser beam has the unique ability to deliver very 
high power per unit area. As the optical energy strikes the 
surface, it is partially absorbed and partially reflected. 
For metal surfaces and metal alloys, the reflectivity is 
typically very high (92% or greater) at room temperature. 
However, for metals interacting with infrared wavelength 
radiation, there is a decrease in reflectivity with 
temperature, with an abrupt drop in reflectivity as the 
material changes from solid to liquid state as shown in 
Table 1. Tgj^ is the melting temperature of the metal. [Ref. 
11.] This loss of reflectivity causes a large increase in 
the amount of energy absorbed, which results in an even 
faster increase in temperature. Provided the amount of 
energy added is greater than the rate of energy diffusion to 
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surrounding cooler areas, the surface will become hotter and 
hotter and change from solid to liquid and finally become 
vaporized. 


Table 1. Absorption Percentage vs. Temp for Metals. (After Ref. 11) 


Metal 

T. 

(K) 

Absorption 

at 300 K 

Absorption 

at T 

xn 

(Solid) 

Absorption 

at T 

XQ 

(Liquid) 

Absorption 

at 

T + 500K 

XQ 

Ag 

1234 

0.4% 

2.0% 

4.3% 

5.4% 

Al 

933 

1.1% 

4.4% 

10.0% 

13.0% 

Au 

1336 

0.6% 

3.5% 

8.3% 

10.3% 

Cu 

1356 

0.5% 

2.8% 

5.8% 

7.3% 


For typical radome materials there are basically two 
mechanisms of destruction. The first is penetration by 
melting and ablating the material. The second is to create 
fractures by thermally stressing the material. 
Additionally, if the missile housed an infrared seeker, the 
radome would naturally be transparent in the infrared. In 
that case, destruction of the radome would not be necessary. 
Instead, the infrared laser would penetrate the radome and 
interact directly with the guidance system and other parts 
behind the radome. [Ref. 12] 
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Some of the physical processes that govern the damage 
caused to the material by the laser include absorption and 
reflection, density of the material, heat capacity, and 
thermal conduction. Furthermore, the characteristics of the 
laser such as irradiation wavelength, power density, peak 
power, average power, and pulse characteristics play an 
important role as well. 

1. Required Power to Destroy Target 

One of the major reasons behind the experiments 
discussed in this chapter was to determine just how much 
power from a short-pulse FEL is needed to destroy a missile 
in the few seconds allowed for an engagement. 

One can estimate the required irradiance by assuming 
that the laser burns through the material by breaking the 
cohesive bonds of individual atoms allowing the removal of 
the atoms. In actuality, the matter would most likely 
disintegrate in "chunks" of atoms. When the laser beam 
strikes the missile, not all the energy is absorbed, and 
even some of the absorbed energy makes its way to cooler 
parts of the material prior to melting. Additionally, some 
of the power may be wasted by heating some of the atoms 
beyond the temperature needed to remove them. 

One of the materials irradiated during the experiments 
was aluminum. Metals such as aluminum could be expected to 
make up the body casing of a cruise missile. The binding 
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energy of aliominiun is approximately 3.4 eV/atom [Ref. 13, 
pg. 74] . If the casing is assumed to be 3 layers thick, 
with each layer approximately 1 cm thick, and the laser spot 
size on the target is 100 cm", the voltime of material to be 
removed is approximately 300 cm\ The atomic spacing for 
aluminum is approximately 2.5 angstroms [Ref. 13. pg. 98], 
or 2.5 X 10'® cm, which corresponds to the total niimber of 
atoms to be removed 


_ 300 cm^ _ 

( 2.5 • 10“*cm)^ / atom 


2 • 10“ atoms. 


( 1 .) 


The total energy required is therefore estimated at 


E = (2 • 10“atoms) (3.4eV / atom) (l. 6 • 10"*’J / ev) 

» lOMJ 

To deliver 10 MJ of energy in 3 to 4 seconds requires 
approximately 3 MW at the target, assxaming all the energy is 
absorbed. Three MW divided over 100 cm" gives a required 
irradiance of 4> = 30 kW/cm". The actual required irradiance 

level would be determined by the amount of reflected energy. 

A second method to estimate the required irradiance is 
to determine the amount of energy needed to bring the 
material to its vaporization temperature. If that the 
energy is delivered at a rate much greater than the heat 
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loss through diffusion, the required energy can be 
determined from 

<D^ = pz (c[t„ - tJ + ah„ + c[% - tJ + AhJ (3 .) 

where is the required flux density, p= mass density, z = 
thickness of material to be burned through, C = Specific 
Heat, T = melting temperature, T = ambient temperature, 

= vaporization temperature, = latent heat of melting, 

and AH^ = latent heat of vaporization [REF. 14, pg. 167]. 
For aluminum the specific values are p = 2700 kg/m , z - 3 
cm, C = 896 J/kg-K, T„ = 855 K, T, = 300 K, T, =2750 K, AH„ = 
4 X 10^ J/kg, and AH^ = 10.8 x 10* J/kg.' For this rough 
estimate the same value for specific heat has been used for 
both the solid and liquid phases, although in actuallity 
these numbers may be different. Using these values in 
equation (3.) gives a required flux density of 

* 1 GJ / m^ (4 .) 

O 

or 

= 100kJ / cm^ (5. ) 

This much energy delivered over a 3 second engagement 
requires an irradiance of approximately <I>o=35 kW/cm , 
corresponding to a total beam power of 3.5 MW in a 100 cm^ 
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spot on the target. This result is consistent with the 3 MW 
estimate developed with the first approximation method. 

A second material irradiated during the experiments was 
Slip-Cast Fused Silica, a furnace tile material. The 
specific parameters are p- 2200 kg/m^, z = 3 cm, C = 920 
J/kg-K, T„ = 1980 K, T„ = 300 K, T^ = 2200 K, AH„ = 1.5 x 10® 
J/kg, AH^ = 2.2 X 10* J/kg. Using these nximbers in equation 

(3.) shows that 35 kJ/cm% corresponding to 10 kW/cm^ over a 
3 to 4 second irradiation in order to burn through the 
material. This is three times smaller than for aluminum 
because the latent heat of vaporization for al\amimam is five 
times that of the silica, and the vaporization temperature 
of aluminum is also higher. The thickness of the Slip-Cast 
Fused Silica used in the experiments was actually only 9 mm, 
about 1/3 of the thickness of the aluminum used in this 
example. Calculating the required energy for the 9 mm 
sample shows that the total energy needed is only 11 kJ/cm^, 
corresponding to a required irradiance of approximately 
Oo = 3 kW/cm^ for a 3 to 4 second engagement. 

These approximations should provide a good estimate of 
the energy that is actually needed to melt through the 
missile. As reported later in this chapter, an irradiance 
level of 10 kW/cm^ was enough to melt through an al\iminum 
sample that was cut to the correct size in order to control 
the thermal diffusion. However, when a sample of Slip-Cast 
Fused Silica was irradiated at the same irradiance level. 
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burn through was only achieved after an extended time of 
almost two minutes. A possible explanation for this will be 
discussed later in this chapter. 

Experiments conducted in the 1970 's and 1980 's in 
conjunction with the MIRACL program indicate that a power 
density in the tens of kilowatts is needed to destroy a 
missile with a dwell time of a few seconds. Experiments 
conducted by the Air Force during the same time period 
determined that a flux of 50 kW/cm" was needed to damage the 
missile structure in less than one second [Ref. 15] . Both 
of these sets of tests seem to validate the order of 
magnitude approximations made in this section. For this 
reason most of the experiments discussed later in the report 
were conducted at an irradiance level of approximately 
10 kW/cm\ 


2. Scaling Laws and Minimal Spot Size 

Since there is no MW-class FEL to perform full-scale 
experiments, scaling is the only way to determine the 
effectiveness of a FEL weapon. As discussed in the previous 
section, a power density of approximately 0 = 10 kW/cm^ over 
a spot size A = 100 cm^ is required to destroy a missile 
with a dwell time of a few seconds. Scaling laws would 
allow predictions of large area damage from small area 
experiments with the confidence that the two were correlated 
correctly. 
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To achieve a power density of <I> = 10 kW/cm^, a 100 W 
FEL must use a spot size of 1 mm^, while a 1 kW laser uses a 
larger spot size of 10 mm\ Scaling of the laser damage 
will only work if the thermal diffusion is independent of 
spot size, which is only true for laser spots larger than 
the characteristic length for thermal diffusion. Schriertpf 
[Ref. 16] calculates the characteristic thermal diffusion 
length to be 


= 27^ 


( 6 .) 


where 



is the thermal diffusivity and 


n K^AT^ 


(7.) 


where t^ is the time required to bring the material from 
ambient temperature to melting temperature, K is the thermal 
conductivity, AT is the temperature change required to 
reach melting temperature, and is the energy flux being 
absorbed by the material. The thermal diffusion length Dj^ 
represents the distance required for temperature T to drop 
to 1/e times its central value. In the semi-infinite 
approximation used by Schriempf, radial heat flow is 
ignored. In order for this to be valid, the spot size must 
be much larger than or the target diameter smaller than 
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If these conditions are not met, heat 


diameter < 

will diffuse away from the laser spot, and the area will not 
be heated effectively. 


3. Pulse Train 

PELS produce short, powerful pulses with a rapid 
repetition rate in the MHz range. This can be corrpared with 
the typical high peak power short pulse lasers that operate 
at very low repetition rates of Hz or kHz. 

When the laser power is delivered in pulses, the peak 
power increases and may cause new effects beyond thermal 
heating. Short pulsed lasers (microseconds and shorter), 
and ultra-short pulsed lasers (picoseconds and shorter) are 
capable or producing impulse damage in addition to the 
thermal damage. The high peak power can cause rapid 
vaporization at the target surface, and subsequent recoil 
from the vapor blowoff can form a strong pressure wave. The 
pressure wave forms a reflected wave at the backside of the 
material, and the superposition of the two waves can 
overstress the material and cause catastrophic failing of 
the material. [Ref. 17] 

Another effect occurs where an ultra-short laser pulse 
delivers energy to a metal at such a fast rate that the 
metal lattice cannot respond, but the electrons can. The 
electrons rapidly increase in temperature so that the 
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difference between electron and lattice temperatures can be 
as much as a few thousand degrees. [Ref. 18] 

Experiments have shown that the damage threshold 
fluence, of a material gets smaller as the pulse length 
gets smaller. Experiments show that there may be as much as 
a factor of 10 advantage when using shorter picosecond 
pulses over the larger 100 nanosecond pulses. [Ref. 19] 

Research has also examined the threshold for ultra- 
short pulses damaging dielectrics. The same trend as was 
observed with metals was observed in the dielectrics. For 
the longer pulses, the damage threshold decreases with 
decreasing pulse length t^ as > down to very short, 
picosecond pulses. [Refs. 20, 21] 

In the experiments conducted at TJTSIAF, the higher 
fluence pulses contained 0.12 J/cm^ of energy, which is 
slightly higher than the ablation threshold for picosecond 
pulses on metals [Ref. 22] . Some ablation of the material 
may have taken place. However, the fluence level was not 
high enough to create the impulse damage and pressure waves 
discussed in the paragraphs above. 

B. HATERIAliS USED IN EXPERIMENTS 

In addition to the aluminum samples tested in the 
experiments, five other materials were also used. The 
materials chosen were either actual missile nose cone 
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materials, or materials that were very similar to missile 
nose cones. The materials are [Refs. 23, 24]: 

1. Phenolic Resin: 

Radome material from a Standard ARM Missile. 

2. Pyroceram: 

A furnace coating type material used in SM-1 
and SM-2 missile radomes. 

3. Slip-Cast Fused Silica: 

Commercially available furnace material. 
Similar to material used as a radome for 
Patriot missile system. 

4. Polyimide Fiberglass: 

A high temperature fiberglass used on 

supersonic missile radomes. 

5. F2 Epoxy: 

Five plies of fiberglass cloth with an epoxy 
binder. This type of material is used in the 
Soviet STYX type missiles. 

1. Thermal Diffusion Calculations 

The first step in preparing for the experiments was to 
calculate the thermal diffusion length for each material. 
The first calculation and test was made for a small piece of 
aluminTom (Al-6061) . This material had density p = 2700 
kg/m^ specific heat C = 896 J/kg-K, thermal conductivity K 
= 180 W/m-K, thermal diffusivity K = 7.44 x 10'^ m^/s. 
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melting temperature T^= 855K, and power density = 10® W/m" 
= 10 kW/cm^ From equation (6.), the thermal diffusion 

length is D = 2 mm, with a time t=0.01 seconds required to 
begin the melting process. In order to melt through an 
aluminum sample, the laser spot must have an area greater 
than 7C(2 mm)^ = 10 mm\ This would make the spot large 
enough so that the heat arriving at the center of the spot 
would not be able to diffuse to the cooler regions outside 
the spot before melting could occur. Alternately, the 
target itself can be made small with the diameter of the 
sample diameter = D so as to prevent the heat from escaping 
to cooler regions of material before heating could take 
place. 

These calculations were experimentally verified with 
samples of Al-6061. A sample with a 1 cm diameter was 
irradiated with a laser spot size of 1 mm^ (much smaller 
than the required 10 mm^) without melting after several, 
minutes. Other samples with a diameter of 2 mm (very close 
to the laser spot size itself) were melted in a few seconds. 

Another target irradiated was Slip-cast Fused Silica 
(SiOj) . A calculation of the thermal diffusion length 
associated with heating the sample to its melting 
temperature of 1980 K was performed using values p = 2200 
kg/m®, C = 920 J/kg-K, K = 1.26 W/m-K, /t = 5 x 10'’ mVs, T„= 
1980K, Oo = 10® W/m® = 10 kW/cm®. Since this material has a 
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low thermal diffusivity, its thermal diffusion length is 
also small, D^=0.02 mm. An insulating material such as 
fused silica, the scaled laser spot can have an area greater 
than 71 ( 0.02 mm)^ = 0.001 mm% which was easily attained with 
the 1 mm^ or larger beams used at TJNAF. 

2. Experimental Results 

The experiments conducted from March 1999 through March 
2000 are reported in great detail by McGinnis, et al. [Ref. 
25] Results from two of the materials are reviewed in this 
section. 


a) Slip-Cast Fused Silica 

One of the more thorough tests was conducted on 
Slip-Cast Fused Silica, SiO^. The front of the sample after 
all irradiations is shown in Figure 1, with the n\americal 
results recorded in Table 2. 
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Figure 1. Slip-Cast Fused Silica 


Table 2. Irradiation Results of Slip-Cast Fused Silica. 


Run 

Number 

Wavelength 

(pm) 

PRF 

(MHz) 

Airflow 

(mph) 

Average 

Irradiance 

(kW/cm^) 

Laser 

Beam 

Radius 

(mm) 

Penetration 

Rate 

(mm/s) 

1 

3.1 

37.4 

83-86 

10 

1.2 

1.32 

2 

3.1 

37.4 

NO 

10 

1.2 

7.5 

3 

4.8 

37.4 

NO 

9 

0.87 

0.2 

4 

4.8 

37.4 

NO 

9 

0.87 

0.081 

5 

4.8 

37.4 

NO 

9 

0.87 

0.2 

6 

4.8 

37.4 

NO 

9 

0.87 

0.125 

7 

4.8 

37.4 

NO 

9 

0.87 

0.081 

8 

4.8 

37.4 

NO 

9 

0.87 

0.35 

9 

4.8 

37.4 

NO 

500 

0.08 

3.0 

10 

4.8 

37.4 

NO 

500 

0.08 

3.0 
































As can be seen in Table 2, there is quite a 
variation in the penetration rates for the different runs. 
As the exposure time was increased, more smoke and melted 
material filled the hole and blocked the path of the laser 
beam causing the penetration rate to decline over time. 
Apparently 9 kW/cm^ is not a high enough irradiance for this 
material to reach vaporization temperature. The heat was 
able to diffuse away quickly enough that the material 
reached a steady state in the liquid form, and would not 
clear for additional penetration. However, runs 9 and 10 
burned through very quickly indicating that 500 kW/cm^ is 
more than adequate to bring the material to vaporization 
temperature. 

Since run four was the only run to burn completely 
through the material at a power density of 9 kW/cm^ it is 
worth a closer look. A digital picture of damage from run 
four was taken through an optical microscope as shown in 
Figure 2. 
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Figure 2. Damage to Slip-Cast Fused Silica Run 4. 

The laser beam diameter narrowed from the front 
surface to the back surface of the hole; its dimensions were 
1.1 mm at the front and approximately 0.7 mm at the back. 
The actual hole (the section of the sample with material had 
been completely removed during lasing), was very close to 
those dimensions ranging from 1.06 mm at the front down to 
0.7 mm at the back. However, damage to the material was 
much more wider, with the melted portion at the surface of 
the sample measuring 5.0 mm in diameter, and the hole at the 
back edge measuring 2.0 mm in diameter. 

Examination of the hole from run four through an 
optical microscope revealed a 1 mm thick layer of melted and 
rehardened SiO^ filling the hole at the back of the sample. 
It was clear from the video and from the rear power meter 
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that burn-through occurred in run four, but melted material 
solidified and resealed the hole at the back of the sample. 

The volume of the hole in run four is estimated 

by, 


V = ;zR^(z)dz, (8.) 

where the radius closely matches the laser beam radius and 
approximately changes linearly as 

R(z) = 0.53 mm - 0.02 z (9.) 

where z is the distance from the front of the sample. Which 
gives a volume of V = 5.6 mm^ or 5.6 x 10’^ cm^ . By a 
similar calculation, the volume of the entire damaged 
region, including the melted and rehardened portion, is 
estimated to be V = 92 mm\ Based on the density of fused 
silica, p = 2.2 grams/cm^, the amount of material removed 
was 0.012 grams, while the amount of material damaged 
(melted) was 0.20 grams. The heat energy deposited during 
run four is given by 


E = Pt = <I)At, (10.) 

where A = 71(0.53 mm)^, 0 = 9 kW/cm^, and the irradiation 

time t = 110 s, so that E = 9.7 kJ deposited. 

The predictions made in section A.l. of this 
chapter indicated that this material could be burned through 
in about 3 seconds at 3 kW/cm% or just over 1 second at the 


23 



9 kW/cm^ used in these experiments. Using those numbers in 
equation (10.) gives an estimation of the energy that would 
need to be deposited to melt through this sample is 

E = Irradiance • Beam Area • time 

= (3 kW / cm^) • ;r ( 5.5 • 10"^cmf • (3 s) (11.) 

= 0.95 kJ 


vice the 9.7 kJ that was actually deposited. 

The amount of energy needed can also be calculated 
using the first method in section A.l. That is by 
calculating the energy needed to break each atomic bond in 
the vol\ime of the SiO, affected. The atomic weight of SiO^ 
is 60 amu, or 1 x 10'"^ grams. Since the density of SiO^ is 
2.2 grams/cm", the number of molecules per cubic centimeter 
is approximately 2 x 10^\ The required energy to break all 
the bonds assuming approximately 4 eV/bond is estimated as 

E = ( 2.2 • 10 “ bonds / cm^) ( 4eV / bond) ( 1.6 • / ev) 

I -L^ • / 

= 12.8 kJ / cm^. 

From the volume calculated from equation (8.), the 
total energy required is predicted to be 

E = (12.8 kJ / cm^)( 5.6 • 10“^cm^) j 

= 0.72 kJ 
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which is very close to the 0-95 kJ predicted in equation 
(11.), but far from the 9.7 kJ delivered during the 
irradiation of run 4. 

The difference in the predicted and actual energy- 
requirements may have come about as a result of the change 
in state of the material. The calculations used in section 
A.l. of this chapter and in the paragraphs above assume that 
the thermal diffusion length is small compared to the beam 
diameter. This assiimption is true for Slip-Cast Fused 
Silica in its normal state, where the thermal diffusion 
length was calculated in section B.l. to be only 0.02 mm. 
However, it is possible that once the material melted, the 
thermal conductivity increased dramatically, increasing the 
thermal diffusion length substantially. This increase in 
thermal conductivity would result in a dramatic increase in 
the rate of heat diffused away from the laser spot, and 
prevent the material from reaching the vaporization 
temperature. This increase seems to be what is taking place 
since even after being irradiated for 110 seconds, the 
material still contains a large volixme of melted and re¬ 
solidified material. 

This postulate can be checked by calculating the 
energy deposited during rxin 9, which was at a much higher 
irradiance and succeeded in melting through the hole in less 
than 3 seconds. The amount of energy deposited for this run 
was 
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( 14 .) 


E = Irradiance • Beam Area • time 

= ( 500 kw / cm^) • ;r( 8,0 • 10”^cm) • ( 3 s) 

= 0.3 kJ, 

again showing an amount of energy much smaller than the == 
10 kJ deposited during burn through for run 4. This 
reinforces the postulate that the energy must arrive at a 
rate high enough to raise the temperature of the material to 
vaporization temperature before the molten material can re- 
solidify. It appears that almost the entire 9.7 kJ of 
energy delivered during run 4 was wasted because it was 
delivered too slowly. 

In a realistic directed energy application, the 
irradiance used in run 9 would not be needed. A much larger 
spot size of 100 cm^ would not allow for significant thermal 
diffusion away from the damaged area. 

The final experiments were conducted in March 
2000. The purpose of the March 2000 experiments was to 
increase the size of the laser beam spot while maintaining 
the same irradiance of approximately 10 kW/cm^. To do this 
the average beam power was increased to 500 Watts. Air flow 
was also added to these experiments to determine if airflow 
across the material would effect the damage process. 

Burn-through of the material occurred for both of 
the runs conducted to the Slip-Cast Fused Silica during this 
experiment. The burn-through time for run 1 was 6.8 sec and 
for mn 2 was 1.2 sec. In rxxn 1, the damage diameter on the 
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front face of the material sample was 4.4 mm, and on the 
back face was 3.1 mm. In run 2 the damage diameter on the 
front face was 5.6 mm and on the back face 1.8 mm. The 
effect of airflow was to reduce the front face damage 
diameter and to increase the burn-through time. 

The volxome of the entire damaged region in run 1 
is estimated by, 


V,T = n ■ RfT{z) • dz , (15 .) 

where z is the depth of the hole. The radius changes 
approximately linearly as 

R,^(z) = 2.2mm - 0.07 • z (16.) 

The volume of the entire damaged region in run 2 
is also estimated by equation (15.), where the radius 
changes approximately linearly as 

R 2 t(z) = 2.8mm - 0.21 • z (17.) 

After doing the above calculations, the volume of 
the entire damage region in run 1 is Vi^=100 mm^ and in run 
2 is V2 t= 105 mm^. The density of the fused silica is p=2.2 
grams/cm\ therefore mass of the entire damaged region in 
run 1 is 0.22 grams and in run 2 is m2^=0.23 grams. 
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The volume of the holes in runs 1 and 2 are 


estimated using radii 


Rjjj(z) = 1.6mm - 0.17 • z 


(18.) 


and 


R 2 h(z) = 0.75mm - 0.04 • z (19.) 

The calculated voliimes of the hole in run 1 is 
V^„=35 mm^ and in run 2 is Vjg=8.9 mm\ The mass of the 
material removed creating a hole in run 1 is mjg=0.08 grams 
and in run 2 is m 2 ^= 0.02 grams. 

These calculations show that the addition of air 
flow had a negligible effect on the overall damaged 
regions, since both have approximately the same volume. 
Air flow was substantial in affecting the amount of material 
removed since the hole voliame increased from V2 h= 8.9 mm^ to 
V 2 h =35 mm^ when the air flow was blowing across the material. 
However, the air flow also decreased the burn through rate. 
So, though the hole was bigger, it took longer to create. 

Comparing the results of the runs 1 and 2 from the 
March 2000 experiments with those of runs 3 through 8 in 
Table 2, it can be seen that shifting the wavelength from 
A = 4.8 nm to X = 3.1 ^tm, and increasing the average power 
from 100 Watts to 500 Watts had a significant effect on the 
penetration rates. Those penetration rates increased 
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dramatically from 0.081 mm/s to 1.32 mm/s with the presence 
of airflow and to 7.5 mm/s without it. 

A possible explanation for the success of the 
March 2000 runs is the increase in beam radius from 0.087cm 
to 0.12cm. This increase in beam area could be enough so 
that even with the larger thermal diffusion length estimated 
in the preceding paragraphs, the new beam area was large 
enough to ensure enough heat was retained to reach 
vaporization. 


b) Polyimide Fiberglass 

Experiments were also conducted using Polyimide 
Fiberglass. The front and back of the sample after all 
experiments were conducted are shown in Figure 3 and Figure 
4. Table 3 shows the data for all the experiments conducted 
on this material over the year long period. Note that runs 
10 through 15 were conducted at 500 W/cm^, while all the 
other runs were conducted at 10 kW/cm\ 

In Figure 4, we see the backside of the sample 

showing that all of the 10 kW / cm^ irradiations completely 
penetrated the sample. (Note: Hole 3 in Figure 3 has 
mistakenly been n\imbered 2 in Figure 4) . As can also be 

seen, none of the 500 W / cm^ irradiance runs (irradiations 
10 to 15 of Figure 3) penetrated the sample. 
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Figure 3. Polyimide Fiberglass (Front View). 
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Figure 4. Polyimide Fiberglass (Back View). 
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Table 3. Polyimide Fiberglass Irradiation Results. 
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During runs 1 through 3, only run 3 actually 
burned through the sample. These three runs have the same 
irradiation parameters with those of runs 7,8, and 9 of 
Table 3 except for the PRF and the wavelength. Comparing 
these runs, it appears that the shorter {A, = 3,lnm) 
wavelength combined with the lower PRF is much more 
effective, resulting in 500% higher penetration depth rate. 
The lower PRF may be more efficient because it has more 
energy per pulse for the same pulse length. 

For the holes that burned through, the presence of 
airflow resulted in larger holes and an increased 
penetration depth rate than those attained without air, 

The larger spot size irradiations conducted with 
500 W average power (runs 16 to 21), caused 3.5 times higher 
penetration depth rate than the ones using 100 W (runs 4 to 
9), even though both were at the same irradiance level, PRF, 
and wavelength. This could be because the larger spot 
allows material to exit the region quicker, or because the 
larger spot better meets the thermal diffusion requirement. 
Investigation with a microscope reveals that there is no 
evidence of melted or rehardened material inside the holes. 


C. CONCLUSIONS 

This chapter describes some of the first measurements 
of laser damage on missile materials from a newly developed 


32 



TJNAF Free Electron Laser. A much more detailed report of 
all the experiments conducted over a year long period has 
also been prepared [Ref. 25]. 

During the experiments conducted on March 12'*’ and March 
23’"'*, 1999, the only parameter changed was the pulse 
repetition frequency, which caused the pulses used in the 
later experiments to have twice the fluence per pulse as the 
pulses on the 12'"’. The bum through rate for the 
experiments of March 23’"* was faster, and it is postulated 
that the increase occurred because the higher fluence 
pulses, containing 0.12 J/cra in each pulse, was slightly 
greater than the ablation threshold for picosecond pulses on 
metal. It is believed that in the experiments of March 12'"’, 
the laser was heating the samples to the melting point, 
whereas during the experiments of March 23’"*, each pulse was 
ablating a small amount of material immediately. [Ref. 22] 

1. Approximation Methods 

Two approximation methods were developed to help 
determine the required energy to be delivered to the target. 
Both methods assiamed that the energy was being delivered at 
a high enough irradiance to overcome any thermal diffusion. 
Results from estimates and earlier experiments tend to 
indicate that an irradiance level of approximately 10 kW/cm^ 
will be required to melt through most missile materials. 
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2. Scaling 

The TJNAF FEL, using several hundred watts of 
continuous average power was used to simulate the damage 
from a MW-class weapon by focusing the beam to a smaller 
spot size to increase irradiance. The goal of these 
experiments was to begin to develop scaling rules that will 
reliably predict the damage of a larger laser without having 
to bare the enormous cost of building the large laser first. 
The experiments showed that the requirement for the laser 
spot to be larger than the thermal diffusion length is a 
valid requirement. The burn through rate increases 
dramatically as the laser beam radius becomes larger than 
the thermal diffusion. Some of the experiments indicate 
that care must be taken when calculating the thermal 
diffusion length, since this value may change as the 
material changes states. 

3. FEL Pulse Format 

The extremely short, sub-picosecond pulse length of the 
FEL beam is a result of the short electron bunches formed as 
they are accelerated in RF cavities. The TJNAF FEL has a 
unique pulse format compared to conventional lasers with a 
rapid sequence of short, powerful pulses. The peak power in 
each pulse is about 100 MW lasting for only about one-half 
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picosecond coming at a rate of 37 MHz. Other studies have 
shown that such ultra-short pulses may give as much as a 
factor of ten advantage in the reduced fluence required to 
produce damage [Ref. 20]. The experiments reported in this 
chapter began to collect data intended to show whether this 
advantage exists. The results indicate that for a fixed 
average power level, a lower pulse repetition frequency with 
the corresponding higher power per pulse provides a faster 
burn through rate. This is probably because the lower PRF 
pulses contained enough fluence per pulse to immediately 
vaporize the material as well as raising it to melting 
temperature. 

As expected there is a strong relationship between the 
penetration depth rate and the laser irradiance. When the 

irradiance changed by a factor of 20 (500 W/ cm^ 10 
kW / cm^) the penetration depth rate changed by almost the 
same factor. 

From earlier studies and the experiments reported in 
this chapter, it appears that there is a possibility of 
decreasing the fluence required to cause damage to a 
material with ultra-short picosecond pulses. This would be 
a most welcome finding for the high energy laser community 
since more power means more fuel, more weight, and a higher 
possibility of thermal blooming of the beam as it propagates 
through the atmosphere. 
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4. Futxire Ea^eriments 

The TJNAF FEL is scheduled for an upgrade to 10 kW of 
power. This increase will allow more flexibility in scaling 
experiments and further tests of scaling itself. Additional 
experiments that include weighing of the sair 5 )les before and 
after each run, new wavelengths, changing wavelengths during 
irradiation, new pulse formats, and other sample materials 
would be beneficial. As experimental procedures are refined 
and the amount of data increases, a more thorough analysis 
of the FEL damage results and comparison to other lasers 
will become possible. Additional experiments can explore 
just how high the laser irradiance needs to be to bring 
various missile materials to vaporization temperature. 

When higher laser power is available in the future, the 
same irradiations should be conducted with bigger spot radii 
in order to coitpare the results and establish a scaling law. 
It is suggested that thicker samples should also be tested 
in order to determine if the penetration rate stays the same 
after the laser beam has penetrated the material several 
millimeters. Samples should be machined to be only slightly 
larger than the laser spot so that thermal diffusion cannot 
occur. Machining the samples to resemble a small nose cone 
shape may also be important, because different angles of 
incidence would be included. 
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Furthennore, a procedure for making the measurements of 
burn through times more accurate should be established, as 
they are very important in determining penetration depth 
rates. 
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III. BASIC PEL THEORY 


The free electron laser (PEL) has three primary 
components: a relativistic electron beam, an undulator, and 
an optical cavity. The energy from the electron beam is 
coupled to the electromagnetic field stored in the optical 
cavity to produce a very high power coherent laser beam. A 
representation of all three primary components of an PEL is 
shown in Pigure 5. 


electron beam drift space 



undulator magnets 


Figure 5. FEL Undulator in a Resonator Cavity 

Pigure 5 also shows a drift space that is included in 
some undulator designs. The drift space allows increased 
gain in the second section of undulator magnets. This 
design, known as a klystron undulator, will be discussed in 
chapter IV. 
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One of the most important qualities of the FEL is- the 
use of a relativistic beam. Since the electrons are moving 
at speeds close to the speed of light, they experience a 
Lorentz contracted undulating magnetic field, which results 
in the emission of electromagnetic pulses at short optical 
or X-ray wavelengths. [Ref. 26] The emitted wavelength can 
be calculated by solving for the resonance condition, and is 
shown later to be approximated by: 



where A = the emitted wavelength, = the length of one 
undulator period, and y = the Lorentz factor of the 
relativistic electron beam. A typical FEL will have = 5 
cm, and y = 100, resulting in an emitted photons with 
wavelength A = 2.5 urn. 

At TJNAF, a 350 keV photo-cathode injector generates 
the electrons. The electrons must be accelerated well above 
that energy to provide for the short optical pulses desired 
for a weapon class FEL. The electrons leave the injector 
and proceed through a RF linear accelerator, where they are 
accelerated to energy levels of 40+ MeV. 


40 


After exiting the accelerators, the electron beam 
enters the magnetic undulator where it interacts with both 
the static magnetic field of the undulator, and the dynamic 
optical field transiting through the undulator at the speed 
of light. The undulator can be helical or linear, with the 
helical undulator producing circularly polarized light and 
the linear undulator producing linear polarized light. The 
electrons are moving at speeds very close to the speed of 
light, and are almost in synchronism with the light wave. 
However, special relativity requires that the finite energy 
electrons travel slower than the light beam. A condition 
called resonance is defined as: the electrons are traveling 
at a speed such that over the distance of one period in the 
undulator, a light wave of wavelength X will overtake the 
electron beam by exactly one wavelength. At resonance, the 
effects of the moving electrons, the optical wave, and the 
static magnetic field of the undulator produce a beat wave. 
This beat wave travels slower than the light wave and is 
called a ponderomotive wave. At resonance the electrons are 
actually moving at the same speed as the ponderomotive wave, 
and the electric field and resulting force the electrons 
experience is constant and determined by the portion of the 
ponderomotive wave that they are traveling with. It is this 
ponderomotive force that causes the electrons to slow down 
or speed up. [Ref. 27] 
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The change in speeds of the electrons depending on 
their phase with respect to the ponderomotive wave causes 
the electrons to bunch. The bunching is critical if the 
electrons are to add coherent radiation to the light beam. 

The undulator is made from either permanent magnets or 
electromagnets arranged to create an alternating magnetic 
field. The polarization can be helical or linear, linear 
being more common. The purpose of the undulator is to 
periodically "wiggle" the electrons as they transverse the 
length of the undulator. This wiggling motion causes the 
electrons to speed up or slow down (absorb or emit 
radiation), depending on which part of the ponderomotive 
wave the electron is in as it is deflected by the undulator. 

As with most lasers, the optical cavity stores the 
light produced by the laser, allowing a fraction to escape 
on each pass of the light wave. The amount of light that 
escapes on each pass is usually determined by one of mirrors 
having a slightly less than perfect reflection coefficient. 
Allowing too much light to escape results in not enough gain 
for the laser to operate, while allowing too little light to 
escape prevents the laser from achieving sufficient output 
power to destroy a target. 
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A. RESONANCE 

Resonance is described above as the condition where 
over the length of one period of the undulator, one 
wavelength of light passes over an electron. Over N such 
periods of the undulator, the light wave should be a 
distance of nA ahead of the reference electron. This 

distance is referred to as the slippage distance, since it 
represents the distance that an electron has slipped behind 
the light wave. From this information the wavelength of the 
light can be calculated. First, the total distance the 
light wave moves ahead of the electrons is given by; 

N/l = (c - vj At 

n; 1 = (c - c/9j At (21.) 

NA = c (l - At 

where N is the number of undulator periods, X is the 
wavelength of the light wave, c = the speed of light, v^ = 
the speed of the electrons along the line of the landulator, 
Pz = v/c, and At = N,^o/v^ = the time for the electrons 

to transverse N undulator periods. Simplifying equation 
(21.) to solve for the wavelength in terms of velocity in 
the longitudinal direction: 
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(22.) 


nA = c (l - 




A = 


1-ygz 


An approximation of the wavelength can be made by using 
a binomial expansion to determine The transverse 

velocity is shown later to be ^, where = eBAo/2mc^ 

Y 

is the dimensionless undulator parameter. Included in 
are e, the electron charge magnitude, B, the undulator root- 
mean—square field strength, A^^ the undulator wavelength, m 
the mass of the electron, and c the speed of light. is a 

dimensionless n;amber that gives a general measure of the 
magnetic field strength. Continuing with the derivation for 
optical wavelength with the velocity of the electrons with 

respect to the velocity of light, fi. 


Y^ = 


1 


1 - 

= 1 


I- PI- PI 




- p. 


(23.) 

(24.) 


+ K^) = 1 - A' 


(25.) 
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= 1 - 


(26.) 


(l + K) 
f 






1 + K , 
2Y^ 


2 ^ 


(27.) 


where the last approximation was made by using the binomial 
expansion since is typically of order unity, and of 
order 10\ Inserting this approximation for back into 
equation (22.), the value for the optical wavelength 
becomes: 


A = 


1 - A 


1 - 


V 

1 + 


A 




1 + K 
2y^ 


/ 

1 + 

1 - 


V 

. V JJ 


(28.) 


The term in square brackets is very small and can be 
dropped from the denominator. However, it is the only 
remaining term in the numerator and must be kept. 
Therefore, the optical wavelength becomes: 
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l+Jl 


(29.) 


X = 


2f 


and since K is usually <1, 




2r 


(30.) 


Equation (29.) shows how the optical wavelength can be 
controlled by three parameters, the undulator wavelength, 
the undulator parameter, which is primarily controlled by 
the strength of the magnetic field, and the energy of the 
electron beam. 

The resonance condition used to calculate the 
approximate wavelength of the light beam is actually not the 
desired state for the electron beam. At resonance just as 
many electrons would be absorbing energy from the optical 
beam and giving energy back to the optical beam. At low 
gain, the optimum gain condition is for the average electron 
phase velocity, v^, with respect to the ponderomotive wave 
to be at = 2.6, which means the electrons are initially 

traveling slightly faster than the ponderomotive wave. This 
condition will cause the electrons to bunch on the positive 
slope of the ponderomotive potential, which will cause them 
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all to slow down (give up energy) together resulting in 
coherent photon emission. 

B. ELECTRON MOTION 

Three fields combine to form the ponderomotive wave, 
producing three forces on the electrons as they move through 
the undulator. Two of the forces are magnetic. They are 
caused by the alternating static magnetic field of the 
undulator, and the moving magnetic field of the optical 
wave. The third force is from the moving electric field of 
the optical wave. For a helical undulator the three fields 
written in CGS units are: 

Bu = (cos(k^z), sin(k(jz), O) (31.) 

Es = Eg (cos((i/’), - sin(^), O) (32.) 

and 

Bs = Eg (sin((^^), cos(^), O) (33.) 

where B^ is the magnetic field of the \indulator, Eg and Bg 
are the electric and magnetic portions of the optical field, 
k is the wave number of the iindulator, and ^=(kz-<yt+^ is 

o 

the phase of the optical wave, where <[) is the initial phase 
at t = 0 and z = 0. 
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The total force acting on the electron is given by the 
Lorentz force equations [Ref. 28, pg. 579]. The equations 
are: 


(is + /? X i) (34.) 

dt me 

and 

^ ^ . is (35.) 

dt me 

Equations (34.) and (35.) are actually four equations. 
Using these four equations plus equation (23.) 


= 


\- PI- PI 


(23.) 


gives a total of five equations and four unknowns. As is 
commonly done, since the motion of interest is in the 


transverse direction, the 


^{rPz) 


dt 


equation will be ignored. 


The transverse components can be written separately as 
components in x and y as: 
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^ = —+ .AxbJ 

dt me 

= —- (E 3 co^w) - ^nPz sin(k„z) - oos{\i/)) (36.) 

dt me 

= — (Es (1 - Pz) cos{\ff) - sin(k^z)) 

dt me 


and 


= — (Esy + A X X bJ 

dt me 

= —t (_ sin(i«/') + cos(k^z) + sin(y/')) (37. ) 

dt me 

^ (- E^(l - ySj sin(^/^) + B„/9, cos(k„z)) 
dt me 


Sinee 1-P^«\, equations (36.) and (37.) ean be 
simplified to: 


= —(- sin(k„z)) (38.) 

dt me 

and 


ArPy) 

dt 


—(BuAcos(k„z)) 
me 


whieh ean be written together as 


(39.) 
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(-sin(k„z),cos(k„z),0) 


(40.) 


dt me 


where the zero has been added to maintain the longitudinal 
dimension coordinate. By inspection equation (40.) can be 
integrated to show that: 


Pr = —— (cos(k„z), sin(k^z),0) 


(41.) 


The second Lorentz equation is: 


dy _ — e 

dt me 


P • Es 


(42.) 


and using equation (32.), this can be written as 


dy 

dt 

dt 


- e 


me 


(A. A) ■ (e. COS yz-'E.^ sin y/.O) 


- e I - K„ 


me 


Eg (cos(k^z) cos yz - sin(kgz) sin yz) 


(43.) 


Now using the trigonometric identity 
cos(a+b)=cos(a)cos(b)-sin(a)sin(b), this equation simplifies 
to 
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(44.) 


^ cos(k^z + 

dt mcy 

^ ^ aosiC + ^) 

dt me}' 


where C = (k + ko)z - 6rt:, the electron phase in the 
ponderomotive wave. 

Next, the time derivative of equation (23.) can be used 
to get a relation between the time derivative of electron 

y , and the time derivative of the fractional 
velocity, , 

y-^{l + kI) = I - Pi 

-2y-^y(\ + K^) = -2P^ A (45.) 

L = A 

y (l + K^) 

With (0 = kc, the electron phase in the ponderomotive 
field and its derivatives are: 


f = (k + k„)z(t) - (Ot 
^ = (k + k„)z - 0) 

^ = (k + k„ )P^c - 0) 

C = (k + kj P, c 


51 



Therefore, 


P 


z 


• • 

(k + k„)c 


(47. ) 


Now, inserting equation (47.) into equation (45.), 


L = 

Y (i + K^) 

L = V* 

Y (l + K^)(k + kjc 


(48.) 


and since k » k^, and kc = 6), this becomes 


• • • 

- = , 

Y (i + K^>y 


(49. ) 


2710 

Now using equation (29.) and replacing a with - and Aq 

0) 


with 


2 ;rc 

co^ 


(0 = 


1 + 


(50.) 


so that equation (49.) becomes: 
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Finally, using equation (44.) for y •. 


• • 

Y = y = cos(^ + <f) 

2^;/ mc)' (52.) 

V* = 2^, ^ cos(^ + 

mcy 


Equation (52.) gives the electron's phase dynamics■ in 
the form of a penduliim equation. It governs the electron's 
phase space motion resulting from the ponderomotive wave. 
Whenever the total phase (^ + ^) of an electron falls 

between ±^/ the cosine term will be positive, resulting 
in an acceleration of that electron. When the phase is 
between ^ and / the acceleration is negative. This 

positive and negative acceleration results in the bunching 
of the electrons on a scale size of the optical wavelength. 
The electrons longitudinal motion is approximately limited 
to ±X . Introducing a the dimensionless variable 


__ ct 
L 

where c is the speed of light, t = time, and L is the length 
of the undulator. The pendulum equation can be rewritten as 
a second order derivative with respect to T as 


C = |a| cos(^ + $) (53.) 

where dr = - dt and lal = . a = |a|e"’' is called the 

c 

dimensionless optical field. From equation (53.) it can be 
seen that the electron motion is a function of the optical 
field. 

C. OPTICAL WAVE EQUATION 

In the absence of the optical field, the electrons 
experience the forces from the undulator field alone. This 
field causes the electrons to \indergo small transverse 
accelerations, resulting in spontaneous transmissions of 
photons in the forward direction. The nuinber of electrons 
in a single pulse is on the order of 10®, and the n\amber of 
photons emitted during the first pass through the undulator 
is enough that a classical field is formed and then 
amplified though stimulated emission. [REF. 29, p. 129] 
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The optical wave equation written in terms of the 


vector potential is: [Ref. 28, pg. 240] 



(54.) 


where Ji is the transverse current density. 

The vector potential for a circularly polarized plane 
wave can be written as a function of z and t only: 

a(z, t) = [cos(^^r), - sin{l^). O] (55.) 

k 


Carrying out the spatial and time derivatives of 
equation (55.), one obtains 


_ 1 

0z^ k dz 


k + 


V 


[cos(<«^),- sin(^),0] + ^ [sin(^«/'), cos(t</'),0] 

dz ) k dz 


+ - ^ fk + [cos(t/),- sin(^«/'),0] + ^ [cos(?/),- sin(</),0] 

k dz * ~ 


E l' 


dz 


dzj k dz^ 

[- sin(<i/'),- cos(^i^),0] 


(56.) 


and 
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(0 d^K dE. 


dt' 

at 


at 


(O 


d(j> 


a^E. 


[codfff)- sin(?i/'),0] + —^ [sin(t«/'), cos(^^/'),0] 


at^ 


- CO 


J 


a^ 

at 

d^d) 

|[cos(v/') - sin(^^^),0] + Eg [cos(</),- sin(v^),0] 


V 

+ Egf-^ - fyl [- sin(^^/’),- cos(^/),0] 


(57.) 


Assxaming that the optical phases and amplitudes vary 
slowly in time and space 


« kEg, « k^, « fyEg, and ^ « Qxj) 

dz oz dt dt 


the second order derivatives are very small and can be 
eliminated. This allows equation (54.) to be written as: 


[■-2 1 a^ 1 

-h \ 

raEg 1 aEg^ 

V 

c^ dt\ 

A\x, tj = 2 

^ dz c dt j 


+ 2 e, 


[cos(^), sin(^^/'), O] 


— + — -^1 [- sin(t«r), cos{i/c), O] = - — J±(x, t) 

az c at ) c 


(58.) 


At this point it is convenient to use the vectors of 
equation (58.) to define two unit vectors: 


€i = [cos(<^), sin(^),0] and €2 = [- sin(^i^), cos(?^/'),0] 
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The wave equation (58.) can now be separated into two 
orthogonal wave equations: 


and 


+ 1 

^dz c dt j 






(59.) 


d(p 

dz 


+ 


1 _ 27r 

c dtj c 


fji 


A ^ 

S2 

J 


(60.) 


Still left to be determined is the transverse current. 
It can be written as the sum of individual electron 
currents: 


J±(x) = -|e|c]^ - ti) (61.) 

i 


where ri is the position of the i'" electron with charge e 

and transverse velocity , and S is the three dimensional 

Dirac delta function. Substituting equation (41.) into 
equation (61.) gives: 


Jiix 


= -|e| c^ —— (cos(k^z), sin{koz), O) - ri) (62 .) 
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which can be projected onto the two unit vectors defined 
above to obtain: 


a. = ?0 (63.) 

iY 

5x^)• a, = Nck^s£ i£fe-±J^^ tf.(; - j,) ,64.) 

i Y 

when the nioinber of undulator periods is very large, N 
»1, then the gain bandwidth curve is very narrow. The 
change in energy of the electrons per pass through the 
undulator is very small and therefore the Lorentz factor can 
be pulled out of the integral. 


. Ji(x) • El = ^ cos(^ + k^z) - ri) (65.) 

Y i 

Jx(x) • f 2 = ^ sin(^ + k„z) <5^(x - ri) (66.) 

Y i 

The sxim of the product of the individual electrons 
times the phase of that electron in the ponderomotive 
potential can be replaced by considering a volume element of 
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electrons of density times the average phase of the 
potential in that voliame element. 



e|cK„/?e 

r 


{cos{y/ + k„z)) 


(67.) 


and 



NcK,p^ 

r 


{sin{y/ + k„z)} 


( 68 .) 


Now, rewriting y/ + k<,z as ^ the two orthogonal 

wave equations (59.) and (60.) become: 


and 


^ 1 

, 3z c 3t , 


2^NKuPe 

Y 


(cos(^ + ip)) 


(69.) 


' r ^ + i. 

*l^3z c 8t^ 


y 


(sin(^ + <p)) 


(70. ) 


Equations (69.) and (70.) can be formed from the real 
and imaginary parts of the following equation: 


a , 1 a V -•* 

9 z c 9t j Y 


(e-«) 


(71. ) 
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Assximing a long electron pulse, there will be no 
spatial dependence in the z direction, and no longitudinal 

d , 

optical inodes are followed. Therefore the — term can be 

dz 

dropped. [Ref. 29, p 132] 


c dt y 



(72.) 


Finally, using the dimensionless variables 



a = 


_ 4;r|e|K„LE3e 


±0 


y me 


and 


8N(e;>lK„L)V^ 

j'^mc^ 


the optical wave equation can be written as 



(73.) 


where a is the derivative of the dimensionless optical 
field, a, with respect to the dimensionless time variable x. 
Equation (73.) shows that with no current, or with zero 
average phase change (no bunching), there is no change in 
the optical field. 
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D. SHORT PULSES 

Recall that at resonance, the slippage distance was 
defined as N/l. Short pulses are defined as pulses whose 

length is less than or equal to the slippage distcince. 

In the previous section, the long pulse assiamptions 
were used, basically stating that because the pulse length 
was larger than the slippage distance, the optical pulse 
continued to see the same electron density throughout the 
transition through the undulator over many passes. For 
short pulses this assiunption cannot be made. 

Two processes are occurring simultaneously during the 
passage through the optical cavity. The first process is 
the electrical pulse slipping backwards with respect to the 
optical pulse. The second is the bunching that is occurring 
on the optical wavelength scale. Due to the slippage of the 
short electron pulses, only the rear section of the optical 
pulse will remain overlapped with the electron pulse by the 
time the electrons have achieved bunching. This results in 
the optical pulse amplification taking place only in the 
rear of the pulse, while the entire pulse experiences the 
losses associated with the cavity. The pendulvim equation 
showed that electron motion (phase change) only takes place 
in the presence of the optical field and is proportional to 
the magnitude of the optical field. 
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a cos(^ + (j)) 


(53.) 




Since the amplification is only taking place on the 
rear portion of the optical pulse, the centroid of the pulse 
grows toward the rear (making the light pulse appear to be 
moving slower than c) with each pass through the undulator. 
Eventually, the slower moving optical pulse is entirely off 
the electron pulse, and no bunching takes place. The 
optical pulse eventually decays due to the cavity losses. 
This effect is shown in the left column of Figure 6. The 
center window of this column shows the optical field, 
a(z,n), growing strong at n = 100, and then decaying again 
around n = 350. 



-1.5 z 1.5-47 V 47-47 V_47 



-1.5 z 1.5-47 V 47 0 n 400 


Figure 6. Simulation of Short Pulse Evolution in the FEL. 
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Figure 6 shows 9 sections arranged in a 3 x 3 grid 
pattern. The lower left comer indicates the position of 
the parabolic electron pulse shape. The electron pulse 
slips back by the slippage distance nX relative to the light 
pulse as both pulses move through the \indulator from t= 0 
to T=l. The middle-left plot shows the evolution of the 
optical field |a(z, n)| as it grows over n = 400 passes through 
the undulator, with the final optical pulse shape shown in 
the top-left window. 

The center graph shows the evolution of the optical 
power spectrum, P(v,n) over n = 400 passes through the 
undulator, where v is the electron velocity or energy. The 
top-center plot shows the final optical power spectrum after 
n = 400 passes. The bottom-middle plot is the weak-field 
gain spectrum, G{v) , showing a maximxam gain occurring at v = 
2 . 6 . 

The lower-right window shows the Power of the optical 

pulse P(n). The middle-right plot shows the evolution of 

* 

the electron spectrum f(v,n) over n = 400 passes, while the 
top-right plot shows the final spectr\am after n = 400 passes 
through the landulator. 

The electron pulse is timed to enter the cavity in 
synchronism with the reflected optical pulse, which has a 
time interval of 2S/c, where S is the one way distance 
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between the mirrors of the cavity. However, as the optical 
pulse starts to move back, the bunching of the electron 
pulse is delayed more and more, causing the amplification 
process to take place even further back on the optical 
pulse. Eventually, the overlap between optical and electron 
pulses is zero, and the optical pulse decays to zero. 

The desynchronism, d, is defined as the displacement 
between the electron and optical pulses at the entrance to 
the optical cavity. It is normalized to the slippage 
distance. d = 0 means the pulses are in perfect 
synchronization. However, since in effect the optical pulse 
is moving back with each pass, the desynchronism needs to be 
adjusted so that the "slowed" optical pulse arrives with the 
electron pulse. This is done by shortening the length of 
the optical cavity by an amo\int AS, resulting in a 

desynchronism value of d = 2AS/N/1. 


E. GAIN 

FELs are generally operated as either amplifiers or 
oscillators. Amplifiers are single pass FELs that reguire 
high gain to develop high power in one pass through the 
vindulator. Oscillators create their power over many passes 
and therefore can operate at low gain. 
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The gain of a FEL is defined as the percentage change 
in the optical field per pass through the oscillator: 


G = 



(74.) 


where ax is the optical field strength at the end of the 
undulator, and a^ is the optical field strength at the 
beginning of the undulator. 

To achieve gain, energy has to be removed from the 
electron beam and added to the optical beam. Therefore a 
basic method to analyze the gain in the optical field energy 


is to determine the change in energy of the electrons. 


A/ 

r 


for each electron, and add these changes together for all 
the electrons. A negative change in average energy of the 
electron beam means a positive change in the energy of the 
optical beam. 

The change in energy of a single electron is 
proportional to its change in phase velocity. The 
instantaneous phase velocity is given in equation (46.) as: 


V 


= ^ = (k + k„)/?,c - (O 


(46.) 


or in units of dimensionless time r where 9t — 8r 

L 
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(75.) 


y = C = - [(k + kj/?,c - co] 
(k + k„)^^ 


V = ^ = L 


V = C = L[(k + kj/?, - k] 


where L = the length of the undulator. 

The change in the phase velocity can be written as: 


Av = A C = L(k + kjAy^, 
Av = A ^ « LkA/3^ 


where the approximation comes from the fact that k » k^. 

To relate Av to Ay use the substitutions L = nAo, k = 

271!X, k^ = 27iIXo, and equation (29.) 


A = 


1 + 
2 /^ 


(29.) 


to obtain: 
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Av 

Av 

Av 


Ni, 

hItA 


^ 27t^ 


\^J 




yX j 


r 


n2;H 


lY 


2 ^ 


1 + K,^ 




(77.) 


Now use equation (48.) in difference rather than 
differential form: 


A / ^ Y^fiz AA 
Y (l + K^) 


to substitute Ay^^ into equation (77.) 


Av 

Av 


n2;H 


2y^ 


1 + k" 


4;zN A)' 


Ar (l + K^) 
Y^ Pz 


Pz Y 


Av 


Am 


Ay 


r 


(78.) 


where the last approximation comes from the fact that 
~ \ . The change in energy of the electron is Ajinc^or 
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* Av 2 

AE, = -me 

^ Am 


(79.) 


The average energy change, AE^, for an electron is: 


AEs = Aj'^^.^iHic 


(Av) 2 
—- me 
Am 


( 80 .) 


where is the siun total change in energy of the 

electron beam from all the electrons. 

Assiaming a monoenergetic beam enters the undulator with 
phase velocity v^, this can be rewritten as: 


- (v) - v„ 2 

AEs = Ay,,,,2mc' = yme 


( 81 .) 


The final term to determine is the average energy term 
for each electron at the end of the undulator, (v) . The 

phase velocity v can be determined by integrating the 
pendulum equation once 


= f C dr = J |a| cos(^ + dr 


( 82 .) 
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In weak fields the changes in the optical field. 


electron phase, and optical phase are small: a ~ 0 and 

^ « 0 . In weak fields the pendulxam equation can be 
rewritten by expanding the phase, in terms of the 
dimensionless optical field term a 


C = ^ +...) = |a| cos(^^'’) + + . . .) (83.) 


dr 


using the zeroth order conditions 


a^'*) = a. 


+ . . . 

f" = fo + v.r + 


(84.) 


To zeroth order the penduliam equation becomes 


00 (®) y ,A\ \ 


(85.) 


To first order in 'a' the equations becomes 


C = cos(^^®^) = a.0 cos(Co + v'or) (86.) 

at 
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which from equation (85.) is 

= ao cos(Co + v^r) (87.) 

Integrating equation (87.) once produces 

^ = v(‘) = ^ [sin(^o + VoT) - sin(^o)] (88.) 

v'o 

Integrating a second time provides the first order 
approximation for 


^(1) = :^ [- cos(^o + ^0^) + cos(^o)] - — [sin(ro)] 

(89.) 

^(1) = _fLj[cos(^o + ^0^) - cos(Co)] + V'ot[sin(Co)l} 

Vo 

Using these equations the second order approximation of 


the pendulum equation is: 


.o(2) 


0.(2) 


. 0 ( 2 ) 


y (f(»> + f") + = |a| cos(f« + f">) 

dr^ 




dr= 


Co + {[cos(^o + Vqt) - cos(^o)] + V(,t[sin(^o)]} + C 


( 2 ) 


''0 
= lal 


cos(C^°^ + 


^0 + ao cos(^o + ^0^) + ^ " ^0 ^ 


(90.) 


The right side can be expanded by the identity 
cos(a+b)=cos(a)cos(b)-sin(a)sin{b) and by the approximations 

cos(^*’^) « 1 and sin(^^‘^) = : 


ao cos(Co + ^ 0 ^) + ^ = ^oh°4Co + ^0^) cos(^<‘>) - sin(Co + ^ 0 ^) sin 

OT 

ao cos{Co + v^r) + ^ = ao[cos(Co + ^ot) - sin(Co + 

OT 

^ (^( 2 )) ^ sin(^o + 

OT 

^ ^ £ 2 . sin(^o + ^o^Xcos(Co + ^ 0 ^) - cos(^o) + sin{Co)} 

OT Vq 

(91.) 


Integrating equation (91.) once gives: 


^(2) ^ 

^0 


1 - Vo^ sin Co cos(^o 


(92.) 
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Equations (89.) and (93.) give higher order 
approximations for the phase and phase velocity of the 
electrons, which can be represented together in phase space 
as shown in the large window of Figure 7. This figure shows 
the phase space position of 20 electrons after transversing 
the distance of the undulator. About half of the electrons 
have bunched together at a phase position n. 



-n/2 ? 3Jt/2 0 T 1 


Figure 7. Simulation of Phase Space Plot of 20 Electrons. 

In Figure 7 the large window shows the electrons 
position, and velocity, v, relative to the ponderomotive 

wave. This example shows only the final phase space 
positions of the electrons. The lower right window shows 
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the evolution in the optical phase, <{»(t), and the upper right 
window shows the gain G(T) . 

Equation (89.) shows that to first order the bunching 
is proportional to the optical field a^, and that some of 
the electrons move forward while some move backwards, 
depending upon their initial phase position, and phase 

velocity, Vq . 

To first order, the final velocity is proportional to 
a,). It is also dependent on initial phase and phase 
velocity, with some of the electrons gaining energy and some 
losing energy. 

To complete equation (81.), the final unknown (v) can 
be approximated as: 


(v) = + . . .) (93.) 

The electrons enter the undulator with a pulse length 
that is many optical wavelengths long. The phase space is 
periodic on the optical wavelength scale, so the average 

change in energy AEs can be determined over one wavelength: 

194.) 

i 0 
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The average for the zeroth order approximation is: 


(v”) = Vo (95.) 

indicating no change in average energy of an electron. 

To first order the average is: 


= Vo + 0 


(96.) 


since the two sine functions in equation (88.) cancel when 
averaged over all values of 

To second order the average is: 



.2 r 


v„ + 0 + 


0 L 


cos(Vo^') - 1 + -^ VoT sin(i'oT) 


(97.) 


which produces a net change in the average energy dependent 
upon the initial phase velocity of the electron. 

Equation (97.) can be placed into equation (81.) to 
give the average phase velocity change in the electron beam. 

The gain in the optical beam was defined as the change 
in energy of the optical beam divided by the initial energy. 
It can be determined on the differential volume element dV 
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level as the change in the energy of the electron beam in 
that volume, divided by the initial optical beam energy in 
the same volume. The energy in a volume element containing 
the optical beam is: 


■'O-beam 


2eI ^ 

—s- dV 

%n 


(98.) 


where E is the Electric field of the beam. 

The average change in energy per electron is given in 
equation (97.). The nvimber of electrons in a volume element 
of light dV is yOFdV, where F is the filling factor, which is 
the ratio of the electron beam area to the optical beam 
area. The gain in the optical beam is therefore: 

(/?Fdv)|rmcH(v) - Vq) / 4 ;Zn] 

2E3dV / 

and using equation (97.) the gain becomes: 


G = 


de. 


e-beaia 


^ j 


G = - 


4;rPgF)'mc a^ 


cos(vor) - I + sin(voT) 


( 100 .) 


Using the dimensionless variable j = 8N(e;zK„L)^PgF / j'^mc^ 
produces 
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G 


( 101 .) 


j[2 - 2 cos(vot) - VpT sin(voT)] 



Equation (101.) shows that in the weak field low 
current approximations used here, the gain is primarily a 
function of initial phase velocity, Vq . Figure 8 below is a 

plot of the gain spectrum for a range of phase velocities 
from - 15 to + 15. Gain G is plotted on the vertical axis, 
while initial phase velocity is plotted on the horizontal 

axis. The computer simulated spectrum shows a peak gain of 
27% when the phase velocity equals 2.6. The spectrum is 
anti-symmetric about Vq = 0 , indicating that the maximum 
loss would occur for a electron beam with a velocity of 

Vo = -2.6. 

0.2691 

0.0 

-0.2691 

Figure 8. Weak Field Gain Spectrum 
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A phase space plot showing the evolutions in phase 
space of 20 electrons with initial phase velocity of 
Vq = 2.6 is shown in Figure 9. The light colored dots 

indicate the phase space positions at the beginning of the 
undulator, and the darker colored dots show the positions at 
the end. As can be seen in the figure, many of the darker 
colored dots are bunched in the phase position of 
approximately ^ — n, and are all moving down in phase 
velocity. This indicates that by the time the electrons 
reach the end of the undulator, they are bunched together on 
the optical wavelength scale, and are losing energy to the 
light beam. 



Gain 0.3 

VA 


0.1 


I I I 1..I.I.1..1.-U-I I.I. I-U4J-J 


-n/2 


3JI/2 0 


Figure 9. Phase Space Evolution for 20 Electrons. 
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F. TAPERED UNDUIiATOR 

The resonant condition was defined in equation (29.) 

as: 


A = 


1 + Ku 




(29.) 


The resonant condition was used to calculate the change 
in energy expression in equation (78.). However, once the 
optical beam amplitude becomes large, the change in energy 
of the electrons is such that they are no longer in 
resonance with the optical beam and the gain decreases. 

One method to maintain resonance and as a result 
maintain gain is to change one of the parameters in equation 
(29.), so that the right side of the equation remains 
constant even with decreasing energy. This can be done by 
decreasing the undulator wavelength, Ao, or by decreasing 
the magnetic field strength resulting in a decreased value 
in the undulator parameter Both methods are referred to 
as tapering the undulator. 

Effectively either one of the changes produces the same 
effect on the electrons as if the electron had been 
accelerated back up to the resonance speed by some 
electrostatic force in the direction of the beam. The 
acceleration of the electrons is given by the pendulum 
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equation, and thus the effect of tapering the undulator is 
included in the modified penduliim equation: 

^ = <5 + |a| cos(^ + <!>) (102.) 

where S is the term representing the new 'acceleration' or 
'torque' that the tapering of the undulator will cause. 

To determine the value of S that will result from 
reducing the undulator wavelength start with equation 

(71.) 


V = i = lXfyi + - k] (71.) 


and calculate the effect that an incremental change in \ 
will have on the phase velocity V . 


V 


= U 


k + 


2n 


LV 





S = 
S = 


* / 2?r 

Av = L (- 1) ^ 

Av = -im — 2 . 


PMo 


(103.) 
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where the last line included L = nAq and ~ 1 • A typical 


value of the tapering is —-— — —0.2 . If N — 100, the value 

for the acceleration is J = -200^(- 0.2) = 40;r . 

The value of 8 achieved by a change in the magnetic 
field strength is determined similarly. Start with equation 
(45.) to relate the undulator parameter to by using a 
binomial expansion: 


r"(l + T^l) = I - Pi 


F 

f 1 + kH 

\ ‘ / 

1 + 




(45.) 


(104.) 


Using this approximation in equation (71.) gives 


V = i = L[(k + - k] 

V = L[ky3, - k + ko/?J 

V = L[k(^2 “!)■*■ ^o^zl 


V = L 


1 - 


1 + K,, 


2y^ 


\ 

f, 1 + kH 

+ kft 

1- 

u 

) 

^ ^ y-l 


II 

- k 

'"l + 

+ ko 



1 J 



(105.) 
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1 + 

where the term - — « 1 was dropped from the second 

parenthesis term, but could not be dropped from the first 
parenthesis term because it was the only remaining piece. 

lL7t 

Continuing by making the substitutions L = , k = —, and 

X 

Z 

using equation (29.) for the value of — the value for the 

A 

acceleration term S is 


= Av = nA 


In I 2 k 


A 2r 


. ^ r 1 \. 

-4;^^ ^kAK„ 


1 + K^2r 

k! \ AK. 


1 L 7 AK., 


1 + K„ 


(106.) 


The amount of tapering that can be applied is not 
unlimited. The general rule is 


|a| > 1^1 < 4^/a > 2n 
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The first recjuirement of |a| > |(^| ensures that some of 
the electrons will be trapped in the ponderomotive potential 
well, known as particle trapping. If a is smaller than 8 
then the particles will all be accelerated in phase space 
out and away from the separatrix. The second requirement 
ensures that the tapering is strong enough to overcome the 
natural untapered deceleration, and the last requirement 
ensures that the system in operating in a strong optical 
field condition where tapering is needed. [Ref. 29, pg. 
157] 


G. EFFICIENCY 

The efficiency Tj of a. Free Electron Laser is defined as 
the fractional loss in energy of the electron beam during 
one passage through the undulator. 

^ (Ar)mc^ ^ (107.) 

where {Ay) is the average energy of an electron at the end 
of the undulator. 

The efficiency of the FEL can be increased by tapering 
the undulator. Assiaming ^ fixed, the effect 

of tapering ± S can be determined starting with equation 
(105.) 
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V = L 


(105.) 



f 1 + ^ 

k - k 


o 

JJ 


Recalling L = nXq this can be written in terms of the 
fixed variables X,Xo, and N as 


V 




2;r 27r 

fi + K] 


2y^ 

K ^ Jj 


(108.) 


Rearranging terms, the equation can be rewritten as a 
solution of 

V = 

= 


For a tapered undulator, + AK„r . Typically, 

AK„ « and V « 2m , so 


I - 


2r^A 

^o(l + Ku) 


1 - 


2m 


(109.) 


83 


Y = 


_ Ao(l + + 2K,,AK,r + (AK,r)^) 


2A 1 - 


2m 


r “ 


2,(l + + 2K,,AK,7) 


1 + 


2m 


2^\ 


1 - 


2m 


A^(l + kI^ + 2K^^AK„r|'l + 


1 + - 

^ 2mj 

. ^ 


2m 


(110.) 


2/1 


where higher order terms involving (AK„T)^and have been 

dropped. Continuing with the approximation, equation (110.) 
can be written as 


, _ >lo(l + Ko) 

^ ~ 2X 


, 2K„„AK„t 

V 

1 . uo u 

1 + n -—r 

1 + - 

(l + Ko). 

iTiN^ 


( 111 .) 


Since the second term in each of the square brackets is 
assxamed small, carrying out the multiplication yields 


+ KL) 

■ 2K,,AK„r V ■ 

1 + ”7-+ 

2A 

(l + KlJ 2mi 


( 112 .) 


At T = 0, y = Yo Slid V' = v^. 


and equation (113.) 


becomes 
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== + <o) \l + J^ 


(113 


As the electrons travel down the undulator, both 
V evolve in T, and away from their initial values j; and . 


y = + Ay 


(114.) 


V = + Av 


(115. 


Using equation (114.), equation (112.) becomes 


f = rl + 2r^Ar + {Ayf 


+ kL) r, 2 K^„Ak T Av 

~ 2k \_ (l + K^J im 2;zn_ 


(116.) 


For N » 1, (AyY « yl, so the smaller term can be 

dropped from equation (116.). Using equation (113.) shows 
that yl can be cancelled from both sides leaving 


y = ^Yo^r = 


_ k^jl + K^J r 2K,,AK„r Av 


(117 
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Multiplying both sides by 



yields 


Ay + K^J 

K Ak r Av 

UO U r 

/ T—r -r 

ro 2r.U 

_(i + K^J 4m j 


Now using equation (106.) to define the taper 5 , 


S = 


f 

-4;2N 

V 


K,: 


1 + K,^ 


>1 


ak, 

Ku 


equation (118.) becomes 


Ay A,{l + K^J 

S + Av' 

Vo ~ 2ylX 

4m 


Using the resonance condition equation (29.) 
y = , equation (119.) becomes 


Ay _ yl [Av - S ' 
Ko rl I - 


For N » 1, y^ ~ Yr’ equation (120.) simplifies to 


( 118 .) 


(106.) 


(119.) 


with 


(120.) 
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( 121 .) 


A}' __ [~ Av - ^ 

L 4;zN 


Equation (107.) can be written as 


(A}') _ - (Av) 

Y Am 


( 122 .) 


where Av = , therefore 


„ ^ ^ - (^f) + K) 

4;zN 


(123.) 


where is the final electron energy, and is the initial 
energy. 

Equation (123.) holds for both positive, ^ > 0, and 
negative, < 0, tapering, and examples of the change in 
efficiency for various taper values will be shown in the 
Chapter V. 
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IV. TJNAF UPGRADE SIMULATIONS 


In order to move from damage experiments with a laser 
spot size of a few millimeters to an actual weapon with a 
laser spot size of about 100 cm% the average power of the 
laser must be increased. The current average power at the 
TJNAF FEL is 1.7 kW, and needs to be increased about 3 
orders of magnitude to several megawatts for weapons 
applications. 

At these higher MW power levels, energy efficiency 
becomes important. To operate more efficiently, the TJNAF 
FEL recirculates the electron beam after giving up energy to 
the light beam. The recirculated electrons are timed to 
arrive back at the entrance to the RF accelerator 180 
degrees out of phase with the alternating electromagnetic 
field. This scheme has two advantages. First, it removes 
energy from the electrons so that they can be sent to a 
final beam dump without risk of producing neutron radiation. 
This requires the beam energy entering the dvimp to be less 
than 10 MeV. Second, the energy removed from the electrons 
can be used to boost the energy of the RF field. A diagram 
of the TJNAF FEL is shown in Figure 10. 
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Recirculation 

adds 

at least 

one 

new constraint 

to 

the 

PEL design. Since 

the 

electrons 

are 

being 

returned 

to 

the 

entrance to the 

accelerator. 

they 

must 

undergo 

turns 


totaling 360 degrees. The turns are accomplished through 
bending magnets. The electron's path through the turns is 
effected by its momentum, and by the magnetic force, which 
is proportional to the magnetic field and the velocity of 
the electron. The magnets are designed so that electrons 
moving through the turn at the correct speed will travel 
through the center of the beam pipe that makes up the turn. 
However, electrons with higher energy will be traveling 
faster than this design speed, and above a certain energy 
level their momentum will carry them into the outside wall 
of the turn, or to different points in the bending magnet 
field, or to different points at the entrance of the 
accelerator. Similarly, electrons with lower energy will 
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turn too quickly, and spread the beam. These effects place 
limits on the spread in energy that can be tolerated for a 
recirculated beam, and will have to be considered in the 
upgrade design. 

As a first step in moving from 1.7 kW power to a 
megawatt class laser, TJNAF has decided to upgrade their 
current design to produce 10 to 20 kW average power. To 
support this goal, simulations were conducted to determine 
gain, power, and induced electron beam energy spread that 
would result from the design parameters chosen by the TJNAF 
personnel. 


A. 20 KW SIMULATIONS 

The first sets of simulations were for a 20 kW design 
and included simulations with and without a dispersive 
klystron undulator section. As shown in Figure 5 of Chapter 
III, an undulator can have a dispersive section that allows 
the electrons to continue transiting the optical cavity 
without any interaction. The effect of this dispersive 
section is to allow the electrons to use their slightly 
altered velocities from the first undulator section to bunch 
electrons prior to entering a second undulator section. 
This results in a significant increase in gain achieved in 
the second undulator section, but causes a reduction in 
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saturated power, as will be discussed in the following 
sections. 

The design goals are to achieve 20kW of average power 
while maintaining a final energy spread of 6% or less to 
ensure that recirculation could take place. In order to 
improve gain in weak optical fields, TJNAF's 20kW design 
included a klystron undulator. The simulations varied the 
dispersive section strength D of this klystron undulator 
[Ref. 9], and varied desynchronism value d between the 
optical and electron pulse in order to compare gain, final 
power, and final electron energy spread. The parameters 
used in the simulations are shown in Table 4. 


Table 4. 20 kW FEL Parameters. 


PARAMETER 


Optical Wavelength A 

1 

q 

Undulator Wavelength 

20 cm 

Peak Undulator Magnetic Field B 

2kG 

Undulator Periods N 

24 

Undulator Length L 

6 m 

Undulator Parameter K 

0.79 

Dispersive Strength D 

0,1,2 or 3 

Electron Energy ymc^ 

200 MeV 

Electron Beam Radius 

0.5 mm 

Electron Pulse Length 

2 ps 

Initial Electron Phase Velocity Spread 

Standard Deviation (Tq 

0.3 

Electron Injection Angle Standard Deviation <7^ 

0.15 

Peak Current I 

200 A 

Average Current 

8 mA 
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2 

Cavity losses (Q = 10) 

10% 

Peak Optical Power within cavity 

40 GW 

Repetition Rate 

74.85 MHz 

Duty Cycle 

4x10'^ 


1. Gain 

Simulations were conducted for values of dispersive 
strength D = 0, 1, 2, and 3, with D=0 representing the 

conventional undulator with no dispersive section. The 
desynchronism between electron pulse and optical pulse was 
also varied from d = 0.0 to 0.30 in steps of 0.02. 

The output from a typical simulation is shown in Figure 
11. The variables used in this particular simulation are 
shown across the top of the figure and are: dimensionless 

current j =2, the electron pulse length = 2, 


angular electron velocity spread cr^ 



0.15 


where A0 is the electron beam angular spread, electron 
phase velocity energy spread = 4;zN ^= 0.3, where 
is the fractional energy spread, desynchronism 


d = ^^^^ = 0.03, resonator cavity quality factor of Q = 10* 

to represent a cavity with practically no losses, undulator 
periods N = 24, and dispersive section strength D = 0. The 
high Q factor is chosen to determine the gain produced from 
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the beam interaction alone, ignoring any loses in the 
resonant cavity. 



-2 z 2 “13 V 13 -13 V 13 



-2 z 2 -13 V 13 0 n 40 


Figure 11. Simulation Results Showing Gain. 


Figure 11 shows 9 sections arranged in a 3 x 3 grid 
pattern. The lower left corner indicates the position of 
the parabolic electron pulse shape. The electron pulse 
slips back by the slippage distance nX. relative to the light 
pulse as both pulses move through the undulator from T= 0 
to T=l. The middle-left plot shows the evolution of the 
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optical field |a(z, n)| as it grows over n = 40 passes through 
the undulator, with the final optical pulse shape shown in 
the top-left window. 

The center graph shows the evolution of the optical 
power spectrum, P(v) over n = 40 passes through the 
undulator. The top-center plot shows the final optical 
power spectrxim after n = 40 passes. The bottom-middle plot 
is the weak-field gain spectrum, G(v), showing a maxim\im 

gain occurring at v = 2.6. 

The lower-right window shows the gain of the optical 
pulse G(n) . The middle-right plot shows the evolution of 
the electron spectrxom f(v) over n = 40 passes, while the 
top-right plot shows the final spectrum after n = 40 passes 
through the undulator. 

The gain results from many simulations at different 
values of desynchronism d is plotted in Figure 12. The 
conventional undulator with D = 0 achieved a maximum gain of 
G = 0.2 at a desynchronism value d = 0.03. The maximum gain 
of G = 0.67 was achieved for the klystron undulator with 
dispersive strength D = 2. This was 3 times higher than the 
gain of the conventional undulator with D = 0. However, the 
higher gain achieved using a klystron undulator comes at the 
expense of a lower power at saturation, as described in the 
next section. 
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Figure 12. Steady State Gains 


2. Power 

A typical output showing the power evolution of the 
optical beam is shown in Figure 13 . All the variables used 
in this simulation were the same as used in the gain 
simulation of Figure 11, except that the simulation has been 
allowed to run to n = 12 00 passes, and the quality factor 
has been reduced to Q = 10. This Q represents the actual 
TJNAF FEL cavity loss of 1/Q = 10%. In Figure 13, the 
optical power as a function of n passes is plotted in the 
lower-right corner, while all the other windows of the plot 
remain the same as described in Figure 11. 
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-5 2 S -47 V 47 0 n 1200 

Figure 13. Simulation Results for Conventional Undulator. 

Figure 14 shows the results of many simulations 
plotting the final average power as a function of dispersive 
strength D and desynchronism value d. Maximum power was 
achieved with the conventional undulator and dispersive 
strength D = 0. The power level achieved was 14 times 

higher than the closest klystron undulator which used 
dispersive strength D = 1. Klystrons with dispersive 

strengths D = 2 and D = 3 resulted in even lower power 
levels. 

The results shown in Figure 14 seem to indicate that 
the conventional undulator is the obvious design to choose 
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in order to achieve maximum power. However, strong fields 
in the conventional undulator cause the trapped-particle 
instability. This causes the electron beam energy spread to 
move outside the 6% design tolerance. 



Figure 14. Final Power vs. Desynchronism 

3. Energy Spread 

Figure 15 shows the final induced electron beam energy 
spread produced by the 4 different undulators for 
desynchronism values d = 0.0 to 0.15. As can be seen in 
this figure, the conventional undulator with D = 0 has a 
peak energy spread of 8.6% occurring at d = 0.007. The 
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energy spread decreases with increasing desynchronism, and 
falls below the desired 6% level at a desynchronism value of 
d = 0.06, and remains below 6% for all larger values of d. 

For the klystron undulator with dispersive strength 
D = 1, the maximum energy spread is 5%, which occurs at 
desynchronism values from d = 0.01 to 0.06. It also shows a 
decrease in the energy spread with an increase in 
desynchronism above d = 0.06. The klystrons with D = 2 and 
D = 3 have an energy spread of about 2% for all values of 
de synchronism. 


Fractional Spread in Energy vs. Desynchronism. 



Figure 15. Energy Spread vs. Desynchronism 
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4. Conclusions of 20kW Simulations 

As expected, the maxiinuin gain was achieved with a 
klystron undulator. Using a klystron undulator with 
dispersive strength D = 2 produced a steady state gain of 
67%. However, this increased gain leads to a tremendous 
reduction in final power. The D = 2 klystron design's final 
power was approximately 15 times smaller than that of the 
conventional undulator. 

The conventional undulator achieved the highest final 
power, which occurred at a desynchronism of d = 0.007. 
However, at this desynchronism, the induced energy spread 
was 8.6%, which is outside the limits imposed by 
recirculation. Increasing the desynchronism to d = 0.03 
reduces the energy spread down to the desired 6%, while 
achieving a final output power that is still 9 times higher 
than the best klystron undulator design. Simulations using 
a conventional undulator with a desynchronism value of d = 
0.03 produced a final power of 19.2 kW, which is close to 
the 20kW objective, and corresponds to a final gain of 19%. 

Results for the 20 kW simulations were presented at the 
20'“ International Free Electron Conference in Williamsburg, 
VA, and were published in Nuclear Instruments and 
Measurements. [Ref. 9] 
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B. 10 KW SIMULATIONS 

Recently, simulations for a 10 kW upgrade design were 
conducted. The design parameters are listed in Table 5. 


Tables. 10 kW FEL Parameters 


PARAMETER 


Optical Wavelength A 

5.8 |im 

Undulator Wavelength 

20 cm 

Undulator Periods N 

24 

Undulator Length L 

4.8 m 

Undulator Parameter K 

1,2,3, or 4 

Dispersive Strength D 

0,0.25,0.5,0.75 or 1 

Electron Energy ymc^ 

150 MeV 

Electron Beam Radius 

0.35 mm 

Electron Pulse Length 

0.5 ps 

Initial Electron Phase Velocity Spread 

Standard Deviation 

0.5 

Peak Current I 

270 A 

Average Current 

10 mA 

Dimensionless Current Density j 

2.9 - 3.7 

Cavity losses 

15% 

Peak Optical Power within cavity 

1.5 MW 

Repetition Rate 

74.85 MHz 

Duty Cycle 

3.75x10'^ 


In addition to conducting simulations for various 
values of dispersive strength and desynchronism, these 
simulations also varied the value of the undulator parameter 
K from 1 to 4 in steps of iinity. Table 6 shows the 
parameters affected when the xindulator parameter is varied: 
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X = the optical wavelength, = the dimensionless pulse 
length, j is the dimensionless current density, and F is the 
filling factor which is the ratio of the area of the 
electron beam to the area of the optical beam averaged over 
the interaction region. 


Table 6. Dimensionless Variables vs. K. 



K=1 

K=2 

K=3 

K=4 

X 

2 ."3 Hm 

5.8 |im 

12 pm 

20 pm 


2.7 

1.1 

0.54 

0.32 

jF 

2.9 

3.7 

3.8 

3.8 

F 

0.06 

0.02 

O.pl 

7xl0‘' 


1. Power 

Simulations were carried our for desynchronism values d 
ranging from 0 to 0.3. The dispersive strengths for these 
simulations were varied from D = 0.0 to D = 1.0 in steps of 
0.25. Figure 16 shows the output for the conventional 
undulator design D = 0, desynchronism d = 0.025, and 
undulator parameter K = 2. 
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Figure 16. Power Simulation for Klystron Undulator 


Figure 17 shows the maximum power achieved for various 
values of K and D = 0 (no klystron) . The maximum power 
was achieved for the K = 2 undulator parameter, with a final 
power of 14.5 kW. This was also the maximiam power achieved 
out of all 360 simulations conducted for the various 
undulator parameters K, dispersive strengths D, and 
desynchronism values d. Figure 17 also shows that the K = 1 
simulations achieved similar final power levels at varying 
desynchronism values. The energy spread for these 
simulations was approximately 5.5%, which was within the 6% 
specification for recirculation. No power is achieved for 
the K = 3 and 4 simulations due to the smaller filling 
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factor F and the shorter dimensionless electron pulse length 


<y^. Note that the actual electron pulse length 
fixed, but the slippage distance nX is larger. 



Figure 17. Final Power, Conventional Undulator 


remains 
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15 


D=0.25 


1 



0.00 0.05 0.10 0.15 0.20 0.25 0.30 

desynchronism, d 

Figure 18. Final Power, Weak Klystron Undulator 

Figure 18 shows the final power after introduction of a 
weak klystron of D = 0.25. The figure shows that the 

klystron has added enough gain to produce power for the K = 
3 simulations, but at a final power level of only 7.5 kW. 
The maximum power achieved using the K = 1 and 2 values has 
dropped from 14.5 kW to approximately 10 kW. The energy 
spread for these simulations was reduced to about 4% because 
of the lower power. 

With stronger klystron values of D = 0.5 and up to D = 
1, all four undulator parameter values produced some output 
power, but less than the desired 10 kW. The maximum average 
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power achieved for the higher D values was for K = 3, D = 
0.50, which produced a peak power of 6.7 kW. 

2. 10 kW simulation Conclusions 

From the 360 simulations conducted for the 10 kW 
design, it appears that the conventional undulator D = 0, 
and the weak klystron undulator with D = 0.25, will meet 
the desired goals. However, both undulator designs have 
marginal areas of operation. For the weak klystron 
undulator of D = 0.25, the energy spread was very good at 
4%, but the final output power was only slightly higher 
than the desired 10 kW. For the D = 0 imdulator (no 

klystron), the output power was much higher at 14.5 kW, but 
the energy spread was about 5.5%, which is close to the 6% 
limit. 

Results for the 10 kW simulations were presented at 
the 22"“ International Free Electron Laser Conference in 
Raleigh, NC, and will be published in Nuclear Instruments 
and Measurements. [Ref. 7] 
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V. 


TAPERING 


It is only near resonance that the coupling between the 
laser and electron beams is significant. As the electrons 
lose energy to the light beam, they move away from 
resonance. In order to maintain or re-establish the 
resonance condition, the undulator properties can be tapered 
so that as the electrons lose energy, the resonant energy 
also decreases. The theory of tapering by either changing 
the undulator wavelength or changing the undulator 
magnetic field strength is discussed in section F of Chapter 
III. 

Without tapering, the electrons begin to bunch early in 
the growing optical field, and eventually move to points in 
phase space where they are no longer losing energy but are 
gaining energy from the light beam. When the optical field 
reaches this level the FEL is said to have reached 
saturation. The tapered undulator increases gain and 
efficiency in strong optical fields and extends the usual 
saturation limit to stronger optical fields by maintaining 
the resonance condition. As electrons lose energy 7mc% the 


resonance condition. 



would indicate that 
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either reducing the undulator wavelength, Xq, or the 

undulator parameter, K, would restore resonance. 

Since the purpose of tapering is to reduce the 
resonance energy to match the reduced energy of the 
electrons, it would seem that there would be no benefit in 
using an inverse taper where the undulator wavelength, 
is increased, or the undulator magnetic field strength, K, 
is increased. However, simulations show that there are 
several possible advantages to using the inverse taper. 

The electron dynamics in a tapered undulator were 
described by equation (102.) of chapter III as 

^ = S + |a| cos(^ + 0 ) (102.) 

where ^ is the electrons phase in the ponderomotive wave, ^ 
is the optical phase, a = |a|e^^ is the complex dimensionless 

optical field strength, and S is the acceleration torque 
term resulting from tapering. Simulations using both 

positive and negative values of S, corresponding to positive 
and negative taper of K or were carried out and the 

results will be discussed in following sections. 
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A. EFFICIENCY AND ENERGY SPREAD 

When designing a high power recirculating free electron 
laser, both efficiency and induced energy spread are 
important. As defined in equation (107.) of Chapter III, 

efficiency rj is 



(107.) 


and is typically highest in strong fields. However, strong 
fields have the adverse effect of inducing a larger electron 
energy spread along the undulator. For our purposes, energy 
spread is defined as the full spread in phase velocity. 


The difficulty for a recirculating FEL operating in 
strong fields arises because some electrons gain energy 
while some lose energy. This makes the spread in energy 
large but does not necessarily make the efficiency large. 

A second difficulty arises when considering high power 
FELs. The small percentage of electrons that may lie in the 
outer tails of the energy distribution cannot be ignored, 
because for high power FELs these tails represent a lot of 
current. When recirculating high current (for high power). 
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this small fraction of the electron beam in distribution can 
cause serious damage to the vacuum beam pipe, as well as 
radiation problems as the electrons collide with the walls 
of the recirculating path. Consider just 0.01% of a lA 
(average current) beam at 100 MeV energy focussed to a 1 mm^ 
spot; the irradiarice is 1 MW/cm\ The damage would occur 
rapidly. Therefore, when discussing energy spread for 
high current recirculating devices, the electrons in the 
tails of the distribution must be included as in equation 
(124.) . 

B. NO TAPER 

In the conventional undulator with weak initial optical 
field a^ < n, the gain spectrum is shown in Figure 19. It 
is symmetric about resonance with a peak gain of 27% for j = 
2, occurring at an initial phase velocity v^ = 2.6. The 
gain spectrxim width is ~ tt , and the gain is zero at 

resonance v^= 0. 


110 



Peak Gain of 
27% at 

Vo = 2.6 

0.2691 

0.0 

-0.2691 

Figure 19. Weak Field Gain Spectrum 

As the optical field grows, the gain spectrum changes 
in shape. Figure 20 shows the spectrum as the optical field 
grows from a^ = small to saturation at a^ = 40. 

ln(l+G) 

a„=40 
gain =1% 


Figure 20. Gain vs. Vo in Increasing Optical Field ao 

The spectrum is still symmetric about resonance, but 
the peak gain has been reduced to approximately 1% at an 
initial phase velocity of = 10. At saturation, the FEL 
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oscillator wavelength will shift to follow the point of 
maximum gain, and in this case will shift from = 2.6 to 

Vg = 10, which corresponds to a shift of about 5% in optical 
wavelength for N = 25. 

In the electron pendulum phase-space determined by 
Equation (102.), saturation in strong fields occurs when 
electrons become trapped in closed orbits as shown in Figure 
21 . 



-n/2 3n/2 0 T 1 

Figure 21. Phase Space for Conventional Undulator. 
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when the penduliom separatrix height, 4*s/a , becomes 
greater than the phase velocity for peak gain in weak 
fields, Vq = 2.6, an increasingly larger fraction of the 

electron beam becomes trapped in the closed orbit region of 
phase space. At saturation, the FEL interaction induces a 
spread in phase velocities which is slightly larger than the 

height of the separatrix, ~ 4-\/a ~ Stt , or about 10% for N 
= 25. This spread is sufficient to effect the recirculation 
of the electron beam. 

C. POSITIVE TAPERING 

When someone speaks of tapering the undulator, they 
usually mean positive tapering: either the undulator period 
is decreased along the undulator, or the magnetic field 
strength is decreased along the xmdulator, or both. With a 
modest taper, the maximum weak field gain is shifted and 

occurs at initial phase velocity » 2.6 - ^ . As can be 

seen by comparing Figure 20 with Figure 22, the weak-field 
gain is reduced for a tapered imdulator. However, strong 
field gain is increased. For a taper of 6n, the weak field 
maximum gain is shifted to = —7 . In strong fields of 
= 40, the gain spectrum has become flattened and 

distorted, and the maximum gain has shifted to ~ 0. In 
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strong fields the gain has increased from G == 1% in the 


untapered undulator, to G = 2% with taper. 


Peak Weak Field Gain Shifts left by 5/2 
and is reduced from the Untapered Gain. 



Figure 22. Gain vs. Vo in Increasing Optical Field ao 
for a Tapered Undulator. 


The separatrix surrounding the closed-orbit paths in 


phase space is given by [Ref. 29, pg. 157] 


vl = - Q + 2a(sin(fJ - sin(C)) (125.) 


where 


Co = In - cos ‘(- 5 /|a|) 


(126.) 
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In the tapered \indulator, the phase paths are modified 
by the acceleration term 5 in the pendulxim equation. For 
S > a., the cos'^ expression has no solution, and there is no 

separatrix surrounding closed orbit paths. 

The change in the separatrix can be calculated starting 
with equation (125.). With a tapered xindulator, the 
separatrix height is reduced, and has been shown by Massey 
[Ref. 30] to be 


h = 4>/a|^^|q|cos"*(-|q|)+ Vl-q^ j 


(127.) 


where q = . 


When the taper is too large, such that 


8>a^, the electrons phase space has no closed orbits. 
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8 >a„ 

Accelerating electrons never 
enter closed orbits. 


Time evolution is shown 
as electrons change from 
light 
to dark. 



Figure 23. Phase Space Showing 5 > ao. 


In sufficiently strong fields, > S, some of the 

electrons can be trapped in the closed orbits of the 
pendulum phase space. This is illustrated in Figure 24. In 
stronger fields, a^ » S, a large fraction of the electrons 
can be trapped near resonance around ^ ~ n , and continue to 
lose energy to the optical field, as is illustrated in 
Figure 25. 
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6»a„ 

Most electrons 
enter closed orbits. 



Time evolution 
shown 

as electrons 
change from ligfbt 
to dark. 

-20 


Gain 


0.1 


0 

0.05 


3n/2 Ox 1 


Figure 25. Phase Space Showing ao » 5. 


As shown in Figure 26 below, electrons starting near 
the phase ^ ~ 0 are accelerated by both the torque S and the 

strong optical field a, while the electrons starting near 
~ n see the acceleration term 6 and a roughly cancel, 

leaving them trapped in closed orbits. Tapering is 
effective because the electrons near the maximum gain point 
of C ~ ^ trapped, while the electrons near the phase 

point for absorption ^ « 0, are taken away from resonance 
and eventually stop interacting. It is this imbalance 
between the gain and loss mechanisms that leads to gain. 
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Torque term 5 and 
optical Held cancel 
together for electrons 

Torque term 5 and starting near Co“ 

optical field a^add 
together for electrons 
starting near 



Figure 26. Separatrix for Positive Tapered Undulator. 


As can be seen in Figure 27, which shows only the final 
position of each electron, the electrons appear to be have 
an energy spread with about half of the electrons trapped 
near resonance v ~ 0, and the other half untrapped at 

V ~ S + 2Va . In a strong field of a = 40 and with a taper 
rate of S = Btt , the induced spread approximately 

Vg == J + 2-N/a = 15 + 2 V 4 O = 30 . Due to this large energy spread, 

the tapered undulator does not appear to be desirable for 
recirculating the electron beam. 
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-ir/2 ^ 3 ;t/2 0 t 1 


Figure 27. Phase Space for Positive Tapered Undulator. 


D. NEGATIVE TAPER 

In the negative or inverse tapered undulator, the 
acceleration term d < 0. Initial investigations and 

simulations into the use of negative taper indicated some 
favorable results could be obtained in the areas of strong 
field gain and energy spread. 

Figure 28 shows the gain spectrum for an FEL with 
inverse taper rate S = -6n over N = 24 periods with a 
current density of j = 2. These parameters are descriptive 
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of the TJNAF FEL with recirculation of the electron beam. 
In weak fields, a < n, the gain spectrum peak is shifted 

above resonance to = 2.6 — ^ = 12, with a peak gain of 

G = 12%. This gain is comparable to the positive taper case 
shown previously in Figure 22 above. Both positive and 
negative taper have the effect of shifting the gain spectrum 
by Vq « 2.6 - and reducing the peak gain available by 

about 50% when 5 = ihK . 

In stronger fields, the gain spectrum distorts and 
flattens, but in a manner that differs from both the non- 
tapered and positive tapered undulators. As shown in Figure 
28, the peak in the gain spectrum tends to stay near the 
same value with only a slight shift from = 12 to 

= 13.5 as the field grows from weak to strong values. 

This means that in a FEL with inverse taper, the expected 
wavelength shifts only slightly as the laser approaches 
saturation. The peak gain in strong fields does decrease, 
but only down to 3%, which provides a significant advantage 
over the untapered undulator, and even betters the positive 
tapered undulator. This result was unexpected since the 
tapered undulator has long been thought to provide the best 
extension of the FEL to strong fields. 
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Strong Field Gain tends to remain near 
weak field gain maximum Vo. Strong Field 
In (l+<3) but to a level significantly 

. higher than the conventional and positive 

1 tapered undulators. 



Figure 28. Gain vs. Vo in Increasing Optical Field ao 
for a Inversely Tapered Undulator. 


The electron evolution in phase space with inverse 
taper is also quite different than the other two undulator 
types. The separatrix shown in Figure 29 is still given by 
equation (124.) above, but with <5 < 0. Figure 29 is shown 

for a taper rate of ^ = -(in, N = 24 periods, and a current 
density of j =2. 
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Electrons start at Vo = 15. 


Electrons cannot 
enter closed seperatrix, 

and must bunch as_ 

they move around 
separatrix to low 
energy level, resulting 
in coherent radiation. 

-7T/2 C 33T/2 0 T 1 

Figure 29. Phase Space for an Inverse Tapered Undulator. 

The electrons started at = 15 , which provides the 

peak gain in strong fields for this taper. All of the 
electrons start in open phase space paths far from resonance 
and cannot enter the closed-orbit region inside the 
separatrix. Following the open-orbit paths, it is clear 
that the electrons must bunch as they lose energy to travel 
around the closed-orbit region. The bunch continues to 
follow open orbit paths down to a final phase-space site 

around V ~ -20 where they have lost considerable energy. 

More importantly, the induced average energy spread in the 
final electron distribution appears to be relatively small. 
Most of the final electrons are located in a spread of 
Vg = 15 around v = -16, with a tail extending back to 
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resonance. This corresponds to only a 4% average energy 
spread. This is significantly reduced compared to the 
tapered undulator spread of 13% and the periodic undulator 
spread of 10%. 

These results are surprising and appear promising for a 
recirculating FEL. For high power recirculating FELs, it 
is the full energy spread that is important. The electrons 
at the extreme edges of the distribution must be considered. 
Figure 30 provides a phase space view of a undulator with a 
taper of 5 = -5n, showing only the final position of the 
electrons. Most of the electrons are bunched around V = - 
20, there are still some electrons at the top of the 
separatrix around v = 8. When considering the full spread 
of electrons, a comparison of Figure 21 and Figure 27 with 
that of Figure 30 shows that the inverse tapered undulator's 
full spread is = 28. This is only a modest improvement 
over the untapered and positive tapered xindulator with a 
full spread of = 30. 
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Figure 30. Phase Space for a Negative Tapered Undulator. 


E. CONCLUSIONS OF THE CONVENTIONAL AND TAPERED UNDULATORS 

The inverse tapered undulator appears to offer 
advantages over both the periodic and tapered undulators. 
The advantages include larger gain in strong fields, smaller 
wavelength shift at saturation, and a smaller induced 
average energy spread. 

The increase in gain over the untapered undulator in 
strong fields appears to be about 3 to 1, while the increase 
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in gain over the negative tapered undulator is about 1.5 to 

1 . 

While both the untapered and positive tapered 
undulators have large shifts in resonant wavelength as the 
optical field grows to saturation, the inversely tapered 
undulator experiences only a minor shift. The result for 
inverse taper is only a slight shift in the lasing 
wavelength as the optical field grows from start-up to 
saturation. 

The induced full width half maximum energy spread is 
greatly reduced for the inverse tapered undulator. However, 
the full width energy spread is only slightly improved by 
using the inverse taper design. 

Initial results from the inverse taper simulations were 
presented at the 21“ International Free Electron Laser 
conference, and published in Nuclear Instruments and 
Measurements. [Ref. 8] Additional results were presented at 
the 22“^ International Free Electron Laser conference, and 
will be published in Ref. 7. 
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LIST OF SYMBOLS 


a Dimensionless optical field 

o 

a Dimensionless derivative of optical field 
a„ Dimensionless optical field at start of \andulator 
A Area 

B , , Magnetic Field in x, y, or z direction 

X y z 

Undulator Magnetic Field 
P Fractional speed (velocity/speed of light) 

c Speed of light in vacuiim 

C Specific Heat 

d desynchronism value 

D Dispersive strength 

Dj^ Thermal Diffusion length 

8 Dimensionless taper 

5^ Three dimensional Dirac delta function 

e Electron charge 

E Energy 

E^ Electric Field 

£ Unit vector symbol 

F Filling factor 

O Irradiance or Flux density 

Required Flux Density 
(j) Initial phase of optical wave 
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G Optical field gain 

Y Lorentz factor 

Y^ Initial Lorentz factor 
Yr Resonance condition Lorentz factor 
AH„ Latent heat of melting 
AH^ Latent heat of vaporization 
T] Efficiency 

j Dimensionless current 
J Current 

k Optical wave vector 

ko Undulator wave vector 
K Thermal conductivity 
Ku Undulator parameter 

K Thermal diffusivity 

L Length 

X Optical wavelength 

Undulator wavelength 
m Mass 

n Number of passes through the undulator 
N Undulator Periods 

V Phase velocity 

Vj Final Phase velocity 
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Initial Phase velocity 
Energy Spread 
P Power 

Q Resonator quality factor 

r^ Distance in the 'i' direction 

R Radius 

R.^ Radius of the hole for the 'i'"”' run 

R.^ Radius of the entire damaged area for the 'i*"' 

p Mass density 

Electron density 
S Resonator cavity length 

Electron phase velocity spread 
Oe Electron angular spread 
Electron pulse length 
T Temperature 

T^ Melting temperature 
Initial temperature 
Vaporization temperature 
Tp Pulse length 

T Dimensionless variable for length 

V velocity 

V Volume 

V.g Volume of the hole for the rxin 

V., Voliime of the entire damaged area for the 'i'"*'' 


run 


run 
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(D 


Radian frequency 


VP Phase of optical wave 

y Thickness of material 

z Depth of hole 

^ Electron Phase 
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